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ABSTRACT 



NASA's Kepler Mission uses transit photometry to determine the frequency 
of earth-size planets in or near the habitable zone of Sun-like stars. The 
mission reached a milestone toward meeting that goal: the discovery of its 
first rocky planet, Kepler-lOb. Two distinct sets of transit events were de- 
tected: 1) a 152 ± 4 ppm dimming lasting 1.811 ± 0.024 hours with ephemeris 

T [BJD] = 2454964.57375lg:S + ^*0.837495i^;S5 days and 2) a 376 ±9 ppm 
dimming lasting 6.86±0.07 hours with ephemeris T [BJD] = 2454971. 67611^;°^^°+ 
A*45.29485lgoQQ7g days. Statistical tests on the photometric and pixel flux time 
series established the viability of the planet candidates triggering ground-based 
follow-up observations. Forty precision Doppler measurements were used to con- 
firm that the short-period transit event is due to a planetary companion. The 
parent star is bright enough for asteroseismic analysis. Photometry was collected 
at 1-minute cadence for > 4 months from which we detected 19 distinct pulsa- 
tion frequencies. Modeling the frequencies resulted in precise knowledge of the 
fundamental stellar properties. Kepler-10 is a relatively old (11.9 ±4.5 Gyr) but 
otherwise Sun-hke Main Sequence star with T^s = 5627±44 K, = 0.895±0.060 
Mq, and — 1.056±0.021 Rq. Physical models simultaneously fit to the transit 
light curves and the precision Doppler measurements yielded tight constraints on 
the properties of Kepler- 10b that speak to its rocky composition: Mp = 4:.56'ti^l 
Me, Rp = 1.4161^:^^^ R®^ and pp = S.Sl^g g cm'^ Kepler-lOb is the smallest 
transiting exoplanet discovered to date. 



Subject headings: planetary systems — stars: individual (Kepler-10, KIC 11904151, 
2MASS 19024305+5014286) — techniques: photometric — techniques: spectroscopic 
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Introduction 



NASA's Kepler Mission, launched in March 2009, uses transit photometry to detect 
and characterize exoplanets with the objective of determining the frequency of earth-size 
planets in the habitable zone. The instrument is a wide field-of-view (115 square degrees) 
photometer comprised of a 0.95-meter effective aperture Schmidt telescope feeding an array 
of 42 CCDs that continuously and simultaneously monitors the brightness of up to 170,000 
stars. A comprehensive discussion of the characteristics and on-orbit performance of the 



instrument and spacecraft is presented in Koch et al. (2010). The statistical properties of 



the stars targeted by Kepler are described by Batalha et al. (2010) 



In January, 2010, the team announced its first 5 planet discoveries (Borucki et al. 



2010a Koch et al. 2010 Dunham et al. 2010 Latham et al. 2010 Jenkins et al. 2010c) 



identified in the first 43 days of data and confirmed by radial velocity(RV) follow-up. One 
of these - Kepler-8b - shows a clear Rossiter-McLaughlin velocity variation which allowed 



for the measurement of the spin-orbit alignment of the system (Jenkins et al. 2010c). The 



"first five" are all short-period giant planets, the smallest being comparable in size to 
Neptune. Collectively, they are similar to the sample of transiting exoplanets that have 
been identified to dat^ the ranks of which currently hover around 100. The median mass 
of the sample is 0.99 Mj with a lO*'' and 90*^^ percentile of 0.24 Mj and 4.1 Mj. The 
median radius is 1.18 -Rj with a lO*'* and 90*^* percentile of 0.81 Rj and 1.5 Rj. The known 
transiting planets are, statistically speaking, Jovian-like in both mass and size and have 
short orbital periods with a median value of 3.5 days and a lO*'' and 90*^* percentile of 1.5 
days and 8.0 days, respectively. 

In June 2010, Kepler released a catalog of 306 stars with planet-like transit signatures 



^http : / / exoplanet . eu/ 
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(Borucki et al. 2010b). Even if the majority turn out to be false positives, the number of 



transiting planets could plausibly more than double from this pool of candidates. Soon, 
we will leave the realm of small-sample statistics and be able to say something meaningful 
about not only the mass and size distribution but also the dynamical and compositional 
nature of exoplanets. Information about composition will fall from those systems for 
which we can derive not only a radius, but also a mass. Dynamical information will fall 
from multiple-transiting planet systems. Five such candidate systems were included in 



the catalog of Borucki et al. (2010b) and described in more detail by Steffen et al. (2010). 



The discovery of the planets orbiting Kepler-9 (Holman et al. 2010) marked the first 



confirmation of a multiple-transiting planet system. Kepler-9 is a G-type star with two 
Saturn-mass transiting planets in a near 2 : 1 orbital resonance. The system is also the first 
to show transit timing variations. Dynamical models of these variations afford us the means 
of determining planetary mass without the need for RV follow-up. 

After removing the transit signals of Kepler-9b and Kepler-9c, a third transit signature 
was identified in the light curve revealing an additional candidate with a 1.6-day period. 



The planetary interpretation of Kepler-9d was validated (Torres et al. 2010) without 



detection of a Doppler signal. Rather, given the star and transit properties, a matrix of 
possible false-positive scenarios was constructed. After eliminating all scenarios which are 
not consistent with the observables, a false-alarm probability was computed that speaks to 
the likelihood that Kepler-9d is consistent with the planet interpretation. In this manner, 
Kepler-9d was validated with high confidence as a super-earth-size planet with radius 
l-64lo;}4 -Re- 
Here, we attempt to define "super-earth" from a radius perspective by noting that 



the 10 M0 upper limit proposed by Valencia, O'Connell, & Sasselov (2006) corresponds 



roughly to 2 for a planet with no water and low Fe/Si ratio (Zeng & Sasselov 2010) 
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Practically speaking, the term super-earth has been loosely used to refer to all planets 
larger than earth and smaller than Neptune. However, this is a broad domain that captures 
not only rocky, dry planets that happen to be larger than Earth, but also ocean planets 
and mini-Neptunes. These finer distinctions will only be possible with measurements of 



planetary properties better than 5% in radius and 10% in mass (Valencia, O'Connell, & 
Sasselov 2006 , 2007 [Fortney, Marley fc Barnes" 2007 Seager et al. 2007 Grasset, Schneider, 
fc Sotin|[2009| ). Pending a mass determination, there is no information as of yet with regards 
to the composition of Kepler-9d. 

The discovery of a short-period, super-earth-size planet is not surprising. Indeed, one 
of the most interesting aspects of the Kepler candidate sample reported by [Borucki et al. 



(2010b) is the fact that the median of the radius distribution is strikingly different than 
that of the known transiting exoplanets. The Kepler candidates have a median radius 
of 0.30 -Rj - smaller than that of Neptune (0.34 - and candidates as small as 1.5 

Should this distribution survive the process of false positive elimination, we will see 
substantial numbers of short-period super-earths. There is precedence already for transiting 



exoplanet discoveries in the super-earth domain. CoRoT-7b (Leger et al. 2009 Queloz et 



al. 


2009; 


Bruntt et al. 


2010 



i?0, 4.8 ± 0.8 M0 planet orbiting a K-type star. The case for CoRoT-7b is complicated by 



activity-induced RV jitter. Independent analysis of the Doppler measurements by Pont, 



Aigrain, & Zucker (2010) reduces the significance of the detection somewhat (yielding 



2.3 ± 1.8 Mq) and puts the 95% confidence interval between and 5 - still, however. 



within the super-earth domain. And while the Queloz et al. (2009) mass and radius point 



to a rocky composition, the lower mass of Pont, Aigrain, & Zucker (2010) marginally favors 
a water/ice composition. 



GJ 1214b (Charbonneau et al. 2009) is another example of a transiting super-earth (as 
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defined from a mass perspective) at 6.55 M^. Orbiting an M-type star with a period of 1.58 
days, GJ1214b has a radius of 2.68 R^. Its density (1.87 g cm~^), consequently, is closer to 
that of water than that of the Earth. The interior structure of GJ 1214b has been modeled 



as an H/He/H20 planet with a rocky core (Nettelmann et al. 2010). 



Here, we report on the discovery of a super-earth-size exoplanet orbiting the G-type 
Main Sequence star, Kepler-10 (KIC 11904151). At Kp= 10.96, the star is bright enough 
for asteroseismic analysis of its fundamental stellar properties using the high precision 
Kepler photometry. The stellar properties are known to an accuracy that allows us to put 
Kepler- 10b sitting squarely in the rocky domain of the mass-radius diagram. The light 
curve shows two distinct sets of transit events: one at 0.837495]^q;qqoqq5 days, referred to as 
KOI-72.01 and the other at 45.29485+°;°°°^^ days, referred to as KOI-72.02 where "KOI" 
denotes a Kepler Object of Interest. 

The Kepler-10 data acquisition, photometry and transit detection are described in 
Section [2] The statistical tests performed on the Kepler photometry to rule out false 
positives are described in Section |3| and the subsequent ground-based observations, 
including precision Doppler measurements, leading to the confirmation of KOI-72.01 are 
described in Section |4j Throughout the first half of this paper, we refer to each event as 
KOI-72.01 and KOI-72.02. However, in the latter half, we begin to discuss KOI-72.01 
in the context of a confirmed planet and call it out accordingly as Kepler- 10b. The 
KOI-72.02 signal, at this time, has not been confirmed by RV. Many of the false-positive 



scenarios have been investigated via BLENDER analysis (Torres et al. 2010) as described in 
Section |5j KOI-72.02 will require analysis beyond the scope of this paper for validation 
at an acceptable confidence level, since eliminating possible astrophysical false positives is 
more difficult given the possibility of an eccentric orbit. The larger orbital separation of 
KOI-72.02 precludes us from assuming that tidal effects will have circularized the orbit. 
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While both the transit duration and RV observations are consistent with a circular orbit, 
the uncertainties are large enough that a much more comprehensive BLENDER analysis 
will be required to validate the outer planet candidate. Thus, we refer to this candidate 
throughout the paper as KOI-72.02. 

From spectroscopy to asteroseismology, the analyses yielding fundamental stellar 
properties are discussed in Section [6] Section [7] contains a description of the light curve 
plus radial velocity modeling that yield the planet properties of Kepler-lOb (and KOI-72.02 
under the planet interpretation). In the case of KOI-72.02 the absence of an RV signal 
implies an upper mass limit. We look for small, systematic deviations in the transit arrival 
times that could indicate the dynamical interaction of multiple planets orbiting Kepler-10. 
The implications of a null detection are discussed in Section |8} Finally, the properties of 
Kepler- 10b are discussed in the context of the theoretical models that speak to the planet's 
composition. Its placement in a mass-radius diagram suggests a dry, rocky composition. 
Moreover, the high density (8.8I2 9 g cm~^) indicates a large iron mass fraction as discussed 
in Section ini 

Kepler's primary objective is to determine the frequency of earth-size planets in the 
habitable zone. The number of planet candidates identified in < 1 year of photometry is 
fast approaching the thousands. Time on the large telescopes required for precision Doppler 
measurements is not only costly, it is insufficient for confirming earth-size planets in the HZ 
of sun-like stars. Instruments are not yet capable of yielding the cm s~^ precision required 
for confirmation of Kepler's most interesting candidates. For the near-term future, the 
team will focus its efforts on quantifying the false-positive rate well enough that it might 
be applied to the collective sample of planet candidates. The observations of Kepler-10 are 
part of that effort. The discovery of Kepler-lOb marks an important milestone for the team: 
Kepler's first rocky planet and the smallest transiting exoplanet discovered to date. 
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2. Kepler Photometry 

The discovery of the planet orbiting Kepler-10 begins with the high-precision Kepler 
photometry. Indeed not one but two periodic transit events were identified in the hght 
curve, producing pipeline statistics that initiated the cascade of verification and follow-up 
efforts leading to confirmation and characterization of Kepler- 10b. 



2.1. Data Acquisition 



The Kepler photometer is a 0.95-m effective aperture, wide field of view Schmidt 
camera in an Earth-trailing orbit. It is designed to yield 20 parts per million relative time 
series precision in 6.5-hours for a 12th magnitude G2 star. The focal plane is comprised of 42 
1024x2200 pixel science CCDs arranged together in 21 roughly square modules covering 115 
square degrees of sky (3.98 arcsec per pixel). Each pair of CCDs forming a module shares a 
common sapphire field-fiattener lens. The coatings deposited onto the field-fiattener lenses 
(and, to a lesser degree, the optics and quantum efficiency of the CCDs) define the effective 
bandpass of the otherwise filter-less photometer, yielding a mean transmission of 52.6% 
between 423 and 897 nm (defining the 5% transmission points). Each CCD requires its own 
transmission function for meaningful interpretation of color-dependent behaviors that affect 
planet characterization (e.g. limb darkening). A description of the instrument is given in 



Van Cleve & Caldwell (2009) and Argabright et al. (2008) while an overview of its in-flight 



performance is presented in Caldwell et al. (2010) and Jenkins et al. (2010b). 



Each CCD is electronically divided into two 1024x1100 output units deflning a 
total of 84 readout channels. While the photometer points at a single fleld continuously 
throughout its heliocentric orbit, it rotates about the optical axis once every ~3 months 
(hereafter referred to as a Quarter) in order to keep the solar panels facing the Sun. 



- 10 - 



Consequently, every star spends each quarter of the year on a different channeL Kepler-10 
(RA= 19^02"^43^05, DEC= +50°14'28".68) falls on channels 36 (module 11, output 4), 80 
(module 23, output 4), 52 (module 15, output 4), and 8 (module 3, output 4) in quarters 
0/1 (spring), 2 (summer), 3 (fall), and 4 (winter), respectively. 

CCDs are read out every ~ 6.5 seconds (6.01982 second integration and 0.51895 
second read time), and every 270 readouts are co-added onboard to form 1765.5 second 
(~ 29.4 minute) integrations (long-cadence). Data for up to 170,000 stars are recorded 
at long-cadence (LC) while data for up to 512 stars are also co-added to a 58.85 second 



1 minute) integration (9 readouts) termed short-cadence (SC) as described in Jenkins 



et al.| ( |2010b[ ); |Gilliland et al.| ( |2010[ ). The LC photometry of Kepler-10 used in the 
analyses reported here was acquired between 02 May 2009 and 09 January 2010 - quarters 
0/1 (spring) up through the first month of quarter 4 (winter). In early January, the 
module containing channel 8 experienced a hardware failure that was not recoverable 
and observations of Kepler-10 were cut short until the subsequent spacecraft roll. Over 
11,000 long-cadence observations are used in this analysis. Short cadence data were also 
collected between 21 July - 19 August 2009 (one month of quarter 2) and between 18 
September 2009 through 09 January 2010 (quarter 3 and the first weeks of quarter 4). The 
short-cadence data were vital in determining the fundamental stellar parameters from an 



asteroseismic analysis (p-mode detection) described in Section 6.2 Approximately 200,000 
SC observations were collected in this time period. Both LC and SC data are used in the 
light curve modeling carried out to characterize the planets. Observations of Kepler-10 
continue in the three quarters each year when the target is not on the failed module. 
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2.2. Light Curves 



Raw flux light curves are extracted by performing an unweighted sum of calibrated 



pixels that have been subjected to cosmic ray removal and background subtraction (Jenkins 



et al. 2010a). The pixels used in the sum are those defining the optimal aperture - the 



set of pixels that optimizes the total signal-to-noise ratio (SNR). The optimal aperture 
is dependent on the local Pixel Response Functioiij^ (PRF), measured on-orbit during 



the commissioning period (Bryson et al. 2010), the distribution of stellar flux on the sky 



near the target (crowding), and differential velocity aberration. A complete discussion of 



Kepler's aperture photometry pipeline (PA) is given in Twicken et al. (2010a). Systematic 



errors, outliers and intra-quarter discontinuities are removed by cotrending against ancillary 
data products in the Pre-Search Data Conditioning (PDC) pipeline module as described 
Twicken et al. (2010b). Figure [l] shows the raw (upper panel) and corrected (lower 



m 



panel) flux time series for Kepler-10. Vertical lines denote the boundaries between quarters. 
Intra-quarter fluxes were normalized by their median flux in order to reduce the magnitude 
of the flux discontinuities between quarters. The largest systematic errors are the long-term 
drifts due to image motion (differential velocity aberration) and the thermal transients after 
safe mode events (e.g. that near day 115). After filtering out transit events, the measured 
relative standard deviation of the PDC-corrected, long-cadence light curve is 62 ppm per 
(LC) cadence. An expected instrument + photon noise is computed for each flux in the 
time series. The mean of the per (29.4-minute) cadence noise estimates reported by the 
pipeline is 36 ppm. Both raw (PA) and corrected (PDC) light curves are available at the 



^The Pixel Response Function (PRF) is a super-resolution representation of the distri- 
bution of starlight over the CCD pixels. It includes not only the effects of the instrumental 
optics, but also intra-pixel sensitivity and pointing jitter. 



-VI- 



Multi-Mission Archive at Space Telescope Science Institute MASlj^ 

With regards to the time stamps associated with each photometric flux, we note that 
Kepler's fundamental coordinate system is UTC for all time tags. Spacecraft times are 
converted to barycentric-corrected Julian Dates at the mid-time of each cadence. 



2.3. Transiting Planet Search 

The Kepler-lOb transits were identified by the Transiting Planet Search (TPS) 
pipeline module that searches through each systematic error-corrected flux time series for 
periodic sequences of negative pulses corresponding to transit signatures. The approach 
is a wavelet-based, adaptive matched filter that characterizes the power spectral density 
of the background process (i.e. anything that is not related to the transit signal itself, 
such as stellar variability, instrumental signatures, differential velocity aberration-induced 
photometric drift , etc) yielding the observed light curve and uses this time-variable power 



spectral density estimate to realize a pre-whitening filter to apply to the light curve ( Jenkins 



et al. 2010d). TPS then convolves a transit waveform, whitened by the same pre-whitening 
filter as the data, with the whitened data to obtain a time series of single event statistics. 
These represent the likelihood that a transit of that duration is present at each time step. 
The single event statistics are then combined into multiple event statistics by folding them 
at trial orbital periods ranging from 0.5 days to as long as one quarter (~90 days). The 
step sizes in period and epoch are chosen to control the minimum correlation coefficient 
between neighboring transit models used in the search so as to maintain a high sensitivity 
to transit sequences in the data. 

KOI-72.01 was identified by TPS in each quarter of data with a multiple event 



http : / / archive . mast . edu/kepler 
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statistic > 15a. The long-period transits of KOI-72.02 were identified by manual 
inspection due to the fact that the current version of TPS does not operate on more 
than one quarter of data at a time making detection of events with periods beyond 
~ 30 days necessarily incomplete. Multi-quarter functionality is slated for the next 
software release in early 2011. The transit depth, duration, period, and epoch are 
derived from physical modeling (see Section [?]) using all of the available data. KOI-72.01 
is characterized as a 152 ± 4 ppm dimming lasting 1.811 ± 0.024 hours with transit 
ephemeris of T [BJD] = 2454964.57375l[];S + ^ * 0.837495;;]:^^ days. KOI-72.02 
is characterized as a 376 ± 9 ppm dimming lasting 6.86 ± 0.07 hours and an ephemeris 
T[BJD]= 2454971.676ll[!:[![J^5| + N * A529485tofoo7l days. 



3. Data Validation 

Astrophysical signals mimicking planet transits are routinely picked up by the pipeline. 
The large majority of such false positives can be identified via statistical tests performed 
on the Kepler data itself - tests that are collectively referred to as Data Validation. Data 
Validation for Kepler's first planet discoveries (Kepler-4b, 5b, 6b, 7b, and 8b) is described in 



Batalha et al. (2010). These are gradually being replaced by pipeline software products such 



as those described in Wu et al. (2010). Only targets passing each of these statistical tests 
are passed on to the Follow-up Observation Program team for further vetting, confirmation, 
and characterization. It is at this stage that stars are assigned a ^^Kepler Object of Interest" 
(KOI) number. Kepler-10 was referred to as KOI- 72 throughout the vetting stages. More 
specifically, the short-period event was referred to as KOI-72.01 while the long-period event 
was referred to as KOI-72.02. We will use these identifiers in the subsequent discussions 
of the analyses that led to confirmation. Here we describe the Data Validation metrics 
- statistics which, taken alone, support the planet interpretation for both KOI-72.01 and 
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KOI-72.02. 



3.1. Binarity Tests 

For each event, the even- numbered transits and odd- numbered transits are modeled 
independently using the techniques described in Section [7| The depth of the phase-folded 
even-numbered transits is compared to that of the odd-numbered transits as described 
in 



Batalha et al. (2010). A statistically significant difference in the transit depths is an 
indication of a diluted or grazing eclipsing binary system. Neither of the transit events 
detected in the light curve of Kepler-10 shows odd-even depth differences outside of 2a 



where a refers to the uncertainty in the transit depths reported in Section 2.3 (6 and 9 ppm 
for KOI-72.01 and KOI-72.02 respectively). 

The modeling allows for the presence of a secondary eclipse (or occultation event) near 
phase=0.5 and reports the significance of such a signal. While its presence does not rule 
out the planetary interpretation, it acts as a flag for further investigation. More speciflcally, 
the flux decrease is translated into a surface temperature assuming a thermally radiating 
disk, and this temperature is compared to the equilibrium temperature of a low albedo (0.1) 
planet at the modeled distance from the parent star. There is a marginal detection (~ 2a) 
of a secondary eclipse associated with KOI-72.01. At just 6 ppm, the flux change is not 
severe enough to rule out the planetary interpretation (i.e. that it is due to an occultation 
of a planetary companion). There is no detection of a secondary eclipse associated with 
KOI-72.02. The binarity tests are consistent with the planet interpretation for both transit 
events in the light curve of Kepler-10. 



- 15 - 



3.2. Photocenter Tests 

To check for false positives due to background eclipsing binaries, we study the behavior 
of flux centroids - the center-of-light distribution in the photometric aperture - and how it 
behaves as a function of time, especially comparing images taken during transits with those 
taken outside of transit. We rely primarily on flux- weighted centroids and modeling the 
expected behavior based on the local distribution of point sources from the Kepler Input 



Catalog and supplemented by high spatial resolution imaging (Section 4.2). Neither Speckle 
nor AO imaging reveal any point sources in the photometric aperture that were not already 
listed in the Kepler Input Catalog. 

The study of centroid behavior is complicated by the fact that the Kepler-10 image 
is saturated in all quarters. In Quarter 2, the star is on the edge of saturation, but it is 
apparent that there is some spilling of charge even in this case. Saturation behavior on the 
Kepler focal plane is known to be conservative (e.g. all charge is captured) but has strong 
pixel-to-pixel variation in details such as saturation level and the fraction of flux that spills 
up and down the CCD columns. 

The high-accuracy technique of fitting the Pixel Response Function (PRE) to the 
difference image formed by subtracting the in-transit image from the out-of-transit image is 
inappropriate for saturated targets, because the saturated pixels do not represent the PRF 
of the star, and such mismatches between the pixel data and the PRF cause significant 
position biases in the PRF fit. As an alternative, we use a modified PRF fit technique, 
described below. 

We compute fiux-weighted centroids by creating out-of-transit and in-transit images 
from detrended, folded pixel time series. Separate images are created for KOI- 72. 01 and 
KOI- 72. 02. For each pixel time series, the de-trending operation has three steps: 1) removal 
of a median filtered time series with a window size equal to the larger of 48 cadences or 
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three times the transit duration; 2) removal of a robust low-order polynomial fit; and 3) the 
application of a Savitzky-Golay filtered time series of order 3 with a width of 10 cadences 
(5 hours). The Savitzky-Golay filter is not applied within 2 cadences of a transit event, so 
the transits are preserved. The resulting pixel time series are folded by the transit period. 
Each pixel in the out-of-transit image is the average of 30 points taken from the folded 
time series outside the transit, 15 points on each side of the transit event. Each pixel in 
the in-transit image is the average of as many points in the transit as possible: three for 
KOI-72.01 and eleven for KOI-72.02. In the former case, though there are fewer points per 
transit, there are significantly more transits to draw from (several hundred for KOI-72.01 
compared to just six for KOI-72.02). Consequently, KOI-72.01 yields smaller uncertainties. 

A fiux-weighted centroid is computed for the out-of-transit image and the in-transit- 
image using all of the pixels in the masl^of Kepler-10. This produces row and column 
centroid offsets AR and AC, and the centroid offset distance D = AE? + AC^. 

Uncertainties of these centroids are estimated via Monte Carlo simulation, where a 
noise realization is injected into 48-cadence smoothed versions of the pixel time series for 
each trial. A total of 2000 trials are performed each for KOI-72.01 and KOI-72.02. The in- 
and out-of-transit images are formed using the same de-trending, folding and averaging as 
the flight data. The measured uncertainties are in the range of a few times 10~^ pixels. 

Table [T] shows the resulting measurements of the KOI-72.01 centroids from quarter 1, 
2, 3, and 4 pixel data, along with the Monte-Carlo-based la uncertainties. The centroids 
are converted into centroid offsets and offset distance with propagated uncertainties. We 

^Each star observed by Kepler has a predefined mask that defines the pixel set that is 
downloaded from the spacecraft. The mask changes from quarter to quarter since the star 
falls on a different CCD channel upon quarterly spacecraft rotation. 
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see that while in quarter 1 there is a > 4cr observed offset, in quarters 2 and 3 the observed 
offset is less than la. Quarter 1 had significant thermal and pointing systematics, while 
quarter 2 had significant pointing drift. These systematics were significantly reduced by 
quarter 3. 

Table [2] shows the measured centroid shifts and uncertainties for KOI-72.02 in quarters 
1 through 4. The KOI-72.02 transit uncertainties in quarter 3 are smaller because only one 
transit was observed in quarters 1, 2, and 4 while two transits were observed in quarter 3. 
Generally, the centroid shifts for KOI-72.02 are larger than KOI-72.01. This larger centroid 
shift may be indicative of the KOI-72.02 transit event being due to a background eclipsing 
binary, or the larger shift may be due to the deeper transits of KOI-72.02 (larger flux 
changes induce larger centroid shifts) combined with the complex behavior of pixels at or 
near saturation. To study this question PRF-fit centroids that ignore saturated pixels are 
computed for the KOI-72.02 in- and out-of-transit images, and the resulting centroid shifts 
are computed. These results are shown in Table |3} This technique indicates no statistically 
detectable centroid motion in quarter 3. The uncertainties were computed via Monte-Carlo 
methods similar to those used to compute the fiux-weighted centroid uncertainties. 

While both the fiux-weighted centroid and PRF fitting methods indicate statistically 
significant centroid motion on KOI-72.02 in quarters 1, 2, and 4, we point out that the 
directions of the centroid offsets are inconsistent from quarter to quarter. This indicates 
that these centroid offsets are not likely to be due to a background object in the sky. 

The observed centroids are compared with the modeled centroids computed using point 



sources catalogued in the Kepler Input Catalog (KIC; Latham et al. 2005) out to 15 pixels 



beyond the mask. As described in Sections 4.2.1 and 4.2.2, no additional point sources 



were identified in Speckle or AO imaging. To generate the modeled out-of-transit image, 
the measured PRF is placed at each star's location on the focal plane, scaled by that star's 
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flux. This provides the contribution of each star to the flux in the mask. For each star Sj in 
the mask, the depth dg^ of a transit is computed that reproduces the observed depth in the 
individual pixels. An in-transit image for each Si is similarly but with Sj's flux suppressed 



by 1 — ds^. These model images are subject to errors in the PRF (Bryson et al. 2010), so 
they will not exactly match the sky. Further, a very simple saturation model is applied 
which spills saturation symmetrically up and down the column. 



The modeled centroids (using the quarter 3 mask definition) are presented in Table |4j 
For both KOl-72.01 and KOl-72.02 the offset expected for a transit associated with 
Kepler-10 (and not a nearby star) is smaller than the la uncertainties in D computed 
from the flux-weghted centroid measurements. Modeled transits on other stars in the mask 
predict centroid shift in excess of lOo", which would be readily observed in our flux-weighted 
centroids. We can therefore rule out all known stars in the mask as responsible for the 
transit signal on both KOI- 72.01 and KOI-72.02. 

There is a radius beyond which any star capable of producing the observed transit 
signal would also induce a centroid shift large enough to be detected in the Kepler data 
when comparing in- and out-of-transit images. We estimate this confusion radius by scaling 
the 3(T centroid offset uncertainty by the observed transit depth as described in Section 



4.1.3 of Wu et al. (2010). For this, we use the quarter 3 uncertainty in D. This radius 
of confusion for KOl-72.01 is 1.17 arcsec, and the radius of confusion for KOI-72.02 is 
0.60 arcsec. The volume subtended by this area on the sky can be used in a BLENDER 
analysis (Section [s]) to assess the probability of encountering an eclipsing binary (capable 
of producing the transit signal) in that volume of the Galaxy. However, in the case of 



Kepler-10 the high spatial resolution imaging described in Sections 4.2.1 and 4.2.2 provides 
tighter constraints on the background star population. 
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4. Follow-Up Observations 

Each of the periodic transit signals identified in the light curve of Kepler-10 passes all of 
the Data Validation tests that might indicate the possibility of a false-positive as described 
in Section [3j The star was passed to the follow-up observing team on July 21, 2009 after 
identification and scrutiny of the short-period transit event (KOI- 72. 01). This initiated a 
series of ground-based observations that began with reconnaissance spectroscopy to confirm 
the stellar parameters in the Kepler Input Catalog and identify any obvious eclipsing binary 



signatures (Section 4.1), continued with high spatial resolution imaging to identify nearby 



stars in the photometric aperture (Section 4.2), and ended with high-resolution, high SNR 



echelle spectroscopy with and without an iodine cell to compute stellar parameters, probe 
magnetic activity, measure line bisectors, and make precision Doppler measurements. The 
follow-up observations do not rule out the planetary interpretation for either of the transit 
signatures. However, they only allow for the confirmation and characterization of the 



short-period candidate, KOI-72.01 as discussed in Section 4.3 



4.1. Reconnaissance Spectroscopy 

Spectroscopic observations on medium-class telescopes are acquired before requesting 
precision Doppler measurements. These "reconnaissance" spectra are used to improve upon 
the photometrically-derived stellar classification from the Kepler Input Catalog (Teg = 5491 
K, logs' = 4.47, Ri, = O.983i?0), identify double-lined spectroscopic binaries, and search for 
indications of radial-velocity variations larger than ~ 1 km s~^ that might suggest a stellar 
companion. The objective is continued false-positive elimination. 

Two reconnaissance spectra were obtained with the Hamilton echelle on the Shane 
3-m telescope at the Lick Observatory. High SNR observations were acquired on two 
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successive nights in August 2009, at heliocentric Juhan dates 2455046.771 and 2455047.758, 
corresponding to phases 0.146 and 0.3251 for KOI- 72. 01, and to phases 0.658 and 0.680 
for KOI-72.02. The spectral order covering about 7.0 nm centered on the Mg b lines 
was correlated against a library of synthetic spectra calculated by John Laird using a 
line list prepared by Jon Morse. The stellar parameters for the template spectrum that 
yielded the highest value for the peak correlation coefficient were Tes = 5750 ± 125 K, 
logg = 4.5 ± 0.25, and vsini = O-O^l km s~^ for an assumed solar metallicity. The errors 
quoted for Tgg and logg are half the spacing of the library grid. The two Lick exposures 
were unusually high SNR compared to typical reconnaissance spectra, with a peak value for 
the correlation coefficient of 0.98. Consequently, we were able to estimate the metallicity 
and to interpolate to finer values of the temperature and gravity, obtaining T^s = 5680 ± 91 
K, log^ = 4.33 ± 0.16, and [Fe/H] = -0.09 ± 0.04. The correlation functions for the 
two Lick observations showed no evidence of a composite spectrum, and the two velocities 
agreed within 0.1 kms~^. Thus there was no suggestion of a stellar companion responsible 
for either system of transit events, and the reconnaissance spectroscopy confirms that 
Kepler-10 is a sun-like, slowly rotating Main Sequence star, and supports the planetary 
interpretation for the transit events. Consequently, the star was scheduled for precision 



Doppler measurements (see Section 4.3). 



4.2. High Spatial-Resolution Imaging 

The more complete our knowledge of stellar fiux sources in the photometric aperture, 
the better we are able to assess the likelihood of a blend scenario in the interpretation of 
the transit event. Much of this knowledge comes from the Kepler Input Catalog (KIC) 
which federates point sources from the USNO-B catalog, the 2MASS catalog, and our own 
pre-launch Stellar Classification Program. Identification of point sources within a 1.5" 



radius requires additional imaging. 



4-2.1. Speckle Imaging 

Speckle imaging of Kepler-10 was obtained on the night of 18 June 2010 UT using 
the two-color speckle camera at the WIYN 3.5-m telescope located on Kitt Peak. The 
speckle camera simultaneously obtained 2000 30 msec EMCCD images in two filters: V 
(5620/400A) and R (6920/400A). These data were reduced and processed to produce a final 
reconstructed speckle image for each filter. Figure |2] shows the reconstructed R band image. 
North is up and East is to the left in the image and the "cross" pattern seen in the image 
is an artifact of the reconstruction process. The details of the two-color EMCCD speckle 



camera are presented in Howell et al. (2010) 



For the speckle data, we determine if a companion star exists within the approximately 
2.5 X 2.5 arcsec box centered on the target and robustly estimate the background limit 
we reach in each summed, reconstructed speckle image. The two-color system allows us 
to believe single fringe detection (finding and modeling identical fringes in both filters) if 
they exist and rule out companions between 0.05 arcsec and 1.5 arcsec from Kepler-10. The 
speckle image was obtained with the WIYN telescope native seeing near 0.7 arcsec, and we 
find no Kepler-10 companion star within the speckle image separation detection limits to a 
magnitude limit of 6 mag in R and 4.5 in mag in V below the brightness of Kepler-10. 



4.2.2. AO Imaging 

Near-infrared adaptive optics imaging of KOI-72 was obtained on the night of 08 
September 2009 UT with the Palomar Hale 200" telescope and the PHARO near-infrared 



camera (Hayward et al. 2001) behind the Palomar adaptive optics system (Troy et al. 
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2000D . PHARO, a 1024 x 1024 HgCdTe infrared array, was utilized in 25.1 mas/pixel 
mode yielding a field of view of 25 arcsec. Observations were performed using a J filter 
(Ao = 1.25/im). The data were collected in a standard 5-point quincunx dither pattern (e.g. 
dice pattern for number five) of 5 arcsec steps interlaced with an off-source (60 arcsec East) 
sky dither pattern. Data were taken at two separate times within the same night - 150 
frames using 1.4-second integration times and 150 frames using 2.8-second integration times 
- for a total on-source integration time of 10 minutes. The individual frames were reduced 
with a custom set of IDL routines written for the PHARO camera and were combined into 
a single final image. The adaptive optics system guided on the primary target itself and 
produced measured Strehl ratios of 0.15 at J with a central core width of FWHM = 0.075 
arcsec. The final coadded image at J is shown in Figure [3j 

No additional sources were detected at J within 6.25 arcsec of the primary target. 
Source detection completeness was accomplished by randomly inserting fake sources of 
various magnitudes in steps of 0.5 mag and at varying distances in steps of 1.0 FWHM from 
the primary target. Identification of sources was performed both automatically with the 
IDL version of DAOPhot and by eye. Magnitude detection limits were set when a source 
was not detected by the automated FIND routine or by eye. Within a distance of 1 — 2 
FWHM, the automated finding routine often fails even though the eye can discern two 
sources. Beyond that distance the two methods agreed well. A summary of the detection 
efficiency as a function of distance from the primary star is given in Table [5] 



4.3. Precise Doppler Measurements of Kepler-10 



We obtained 40 high resolution spectra of Kepler-10 between 2009 Aug 31 and 2010 
Aug 06 using the HIRES spectrometer on the Keck I 10-m telescope (Vogt et al. 1994). 
We used the same configuration of HIRES that is normally used for precise Doppler work 
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of nearby FGK stars (Marcy et al. 2008) which yields a Doppler precision of 1.0-1.5 
ms~^ depending on spectral type and rotational vsini. The HIRES fiber-feed was not 
used for these observations. The standard iodine cell was placed in the telescope beam 
to superimpose iodine lines directly on the stellar spectrum. As the iodine lines and 
stellar lines are carried by exactly the same photons hitting the same optics, both sets of 
lines share precisely the same instrumental profile and wavelength scale. The iodine lines 
represent exactly the same spectrometer optics as the stellar spectrum, with no difference. 
We fit the composite spectrum of iodine and stellar lines simultaneously in each 100 pixel 
segment of spectrum, yielding a Doppler shift that automatically includes the instantaneous 
wavelength scale and instrumental profile. This dual-fitting limits the long-term and 
short-term systematic Doppler errors at ~1.0 ms^^. 

Most observations were made with the "C2 decker" entrance aperture which projects 
to 0'.'87 X 14'.'0 on the sky, giving a resolving power of about 60,000 at 5500 A and enabling 
sky subtraction (typical seeing is 0'.'6 — 1'.'2). A few observations were made with the B5 
decker O'.'S? x 3'.'0 that does not permit sky subtraction. It is possible that a few of those 
observations suffered from minor moonlight contamination. The average exposure was 30 
minutes, with some as short as 15 min and others as long as 45 min, depending on seeing 
and clouds. 

The raw CCD images were reduced by subtracting an average bias, subtracting the 
sky counts at each wavelength just above and below the stellar spectrum, flat-fielding 
the spectrum with a 48-exposure sum from a quartz lamp, and extracting the spectrum 
with a width that includes 99.99% of the spectrum. Cosmic rays were removed from the 
raw image first. The pixels typically contained approximately 20000 photons giving a 
Poisson-limited signal-to-noise ratio of 140. We performed the Doppler analysis with the 



algorithm of Johnson et al. (2009). The internal Doppler errors (the weighted uncertainty 
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in the mean of 400 spectral segments) are typically 1.5-2.0 ms~^. The resulting velocities 
are given in Table |6] and shown in Figure |4] as a function of time. The error bars include the 
internal Doppler errors and an assumed jitter of 2 ms~^ (see below), added in quadrature. 
The center of mass velocity relative to the solar system barycenter (Gamma Velocity) for 
Kepler-10 is -98.93 ±0.02 km s"^ (Table (t)). This is an unusually large RV, indicative of old 
disk or even halo membership. The low metallicity (Section 6.1[ ), [Fe/ H] = —0.15 ± 0.04, 



magnetic activity (Section 6.1), and asteroseismic age, 11.9 ± 4.5 Gyr (Section 6.2), also 



suggest old disk or halo membership. 

Nearby stars with Tes near 5600 K and logg near 4.35 such as KOI-72 have been 
previously surveyed for precise Doppler work, revealing a noise-like "jitter" of ~2.5 
ms~^ caused by surface effects including turbulence, spots on the rotating star, acoustic 
oscillations, and atmospheric flows associated with magnetic flux tubes. While each effect 
has its own time scale, it is practical to account for jitter by simply adding it in quadrature 
to the internal errors to yield an estimate of the total uncertainty in the star's velocity. 
We have included a jitter of 2.0 ms~^ in the model fit to all of the data, photometric and 
velocities. 



A periodogram (Figure |5]) of the velocities exhibits a tall peak at a period of 0.837 d, 
in agreement with the photometric period of the KOI- 72. 01. The coincidence between the 
transit and RV periods to three significant digits suggests that the RV period is physically 
related to the transits, as expected if the RV periodicity stems from the reflex motion of the 
star in response to the gravitational influence of the planet. The periodogram also shows 
a peak at a period near 1.2 d, which is the alias resulting from the nightly observational 
cadence. Similarly, there is another peak near a period of 5 days (off the figure) that is the 
alias caused by the beating of the 0.827 d period with the 1.0 d cadence of observations. 
There is no indication of power at the period of KOI- 72. 02. 
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The velocities phased to the photometric period of KOI- 72. 01 (Figure |6| show a clear, 
continuous, and nearly sinusoidal variation consistent with a nearly circular orbit of a 
planetary companion. The lack of any discontinuities in the phased velocity plot argues 
against a background eclipsing binary star as the explanation. Such a binary with a period 
of 0.83 d would have orbital semi-amplitudes of hundreds of kilometers per sec, so large 
that the spectral lines would completely separate from each other, and separate from the 
lines of the main star. Such breaks in the spectral-line blends would cause discontinuities in 
the velocity variation, which is not seen here. Thus, the chance that the 0.83 d periodicity 
exhibited independently in both the photometry and velocities might be caused by an 
eclipsing binary seems quite remote. 

Precision Doppler measurements are used to constrain the mass of KOI-72.01 (Kepler- 
10b) as discussed in Section [7j The absence of a Doppler signal for KOI-72.02 is used to 
compute an upper limit to the mass of this candidate under the planet interpretation. 



4.4. Bisector Analysis 



From the Keck spectra, we computed a mean line profile and the corresponding mean 
line bisector. Time-varying line asymmetries are tracked by measuring the bisector spans 
- the velocity difference between the top and bottom of the mean line bisector - for each 



spectrum (Torres et al. 2005). When RV variations are the result of a blended spectrum 



between a star and an eclipsing binary, we expect the bisectors to reveal a phase-modulated 



line asymmetry (Queloz et al. 2001 Mandushev et al. 2005). In the case of Kepler-lOb 



there is no evidence for a correlation between the bisector spans and the RVs which would 
otherwise argue against the planetary interpretation (see Figure [T]), and similarly for the 
45-day signal of KOI-72.02. However, we note that the uncertainties in the bisector span 
measurements are quite large so that the RMS variation of the bisector spans (10.5 ms^^) 
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exceeds the semi-amplitude of the RV variation (3.315 o ^). Therefore, we do not 
consider the bisector span measurements to be discriminating in this case. 



5. BLENDER analysis of the Kepler light curve 

In this section we examine the possibility that the transit signals seen in the Kepler 
photometry of Kepler-10 are the result of contamination of the light of the target by an 
eclipsing binary along the same line of sight ('blend'). We consider as potential false 
positives physically associated hierarchical triple systems as well as chance alignments 
(eclipsing binary in the background or foreground). We make use of the technique referred 



to as BLENDER, described recently by Torres et al. (2010), which we apply separately to 



each of the signals in Kepler-10 since each could be due in principle to a separate blend. 
Briefly, this technique compares in a sense the observed light curve to a synthetic light 
curve resulting from brightness variations of an eclipsing binary being attenuated by the 
(typically) brighter star Kepler-10. The parameters of the eclipsing binary are varied over 
wide ranges to find all viable blend scenarios producing a good match to the observations. 
The properties of each component of the binary (referred to here as the 'secondary' and 



'tertiary') are taken from model isochrones (Marigo et al. 2008), and those of the main 



star (the 'primary') are constrained by the spectroscopic analysis described earlier. For the 
technical details of BLENDER we refer the reader to the previously cited work, as well as 



Torres et al. (2004). 



5.1. KOI-72.01 (Kepler-lOb) Signal 



Given the short period of this signal, we may assume that tidal forces have circularized 
the orbits of any potential eclipsing binary contaminants (Mazeh 2008). We considered 
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first the case of a hierarchical triple system. Simulations with BLENDER clearly indicate that 
such systems in which the eclipsing binary is composed of two stars provide poor fits to the 
Kepler light curve. We thus rule out this type of blend scenario. If the eclipsing binary 
is composed of a planet (i.e., a smaller tertiary) transiting a star, rather than two stars 
echpsing each other, then it is possible to reproduce the measured light curve, but only if 
the secondary has very nearly the same brightness as the target star itself. In that case the 
resulting size of the tertiary is ^/2 larger than in a model of a single star transited by a 
planet. However, such a bright contaminant would have been evident in our spectroscopy 
as a second set of lines, and this case is therefore also excluded. 

We next examined the background eclipsing binary scenario, allowing the relative 
distance between the binary and the main star to vary over a wide range. We accounted for 



absorption from dust along the line of sight as described by Torres et al. (2010), adopting a 
representative coefficient of differential extinction of = 0.5 mag kpc~^. Interestingly, we 
found that no combination of relative distance and stellar properties for the eclipsing binary 
(composed in this case of two stars) gives an acceptable fit to the light curve. The reason 
is that all such blend configurations lead to out-of-eclipse brightness changes (ellipsoidal 
variations) with an amplitude so large as to be ruled out by the data. Thus, background 
blends of this kind can be confidently ruled out. This result is significant, because it reduces 
the overall likelihood of blends for KOI-72.01 considerably, as we describe below. If we allow 
the tertiary in the eclipsing pair to be a planet instead of a star (i.e., an object of smaller 
radius), then we do find a range of blend scenarios that lead to acceptable fits to the light 
curve, which cannot be ruled out a priori. This is illustrated in Figure [8} in which we show 
contours of equal goodness of fit of the light curve compared to a standard transit model fit. 
BLENDER indicates that these false positives can be up to about 5 magnitudes fainter than 
the primary in the Kepler band, and that in all cases the secondary star is close in spectral 
type (and mass, or color) to the primary, or slightly earlier. Of these blends, we can further 
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rule out those with secondaries having ^Kp < 2, which are bright enough that they would 
have been detected spectroscopically. This implicitly places a lower bound also on the size 
of the tertiaries (~3.8-R0), as tertiaries smaller than this limit only give good matches to 
the light curve if the contaminating star-planet pair is not too far behind the primary, 
and is therefore relatively bright]^ Those cases would be ruled out spectroscopically, as 
mentioned before. The remaining blends shown in the figure above the line with AKp = 2 
must be addressed statistically. 



For this we followed closely the methodology applied by Torres et al. (2010) for the 



case of Kepler-9d. We computed the mean density of stars (i.e., background contaminants) 
in the appropriate mass range based on Figure |8] in half-magnitude bins, using the Besangon 



Galactic structure models of Robin et al. (2003), and we calculated the fraction of these 



stars that would remain undetected after our high-resolution imaging observations described 
earlier (the constraints from centroid motion analysis are less stringent; see Figure [9]). Some 
of these stars might be orbited by (transiting) planets, constituting potential blends. To 
estimate how many of these cases one would expect, we adopted the same frequencies of 



transiting Jupiters and transiting Neptunes as in Torres et al. (2010), based on the results 



of Borucki et al. (2010b), adjusted in the case of the Neptune-size planets to account 



for the lower limit of ~3.8-R0 allowed by BLENDER for the tertiaries. The outcome of 
these calculations is presented in Table [8| The total frequency of false positives (blend 
frequency, BF) we expect to find a priori for KOI-72.01 is BF = 1.4 x 10~^. Translating 



this into a probability statement for the planet likelihood is difficult, as argued by Torres 



et al. (2010), because it requires knowledge of the rate of occurrence of super-earth-size 



^This lower limit of 3.8 also excludes white dwarfs as possible tertiaries. Addition- 
ally, such massive objects in a tight 0.84-day orbit would lead to very significant ellipsoidal 
variation due to tidal distortions induced on the primary star, which are not seen. 
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planets, a quantity that is not yet in hand. Nevertheless, following those authors we may 
express the false alarm rate for a random candidate star in the Kepler field generally as 
FAR = Npp/{Np-p + Np), in which A^fp is the number of false positives and Np is the 
unknown number of planets in the sample. The number of false positives may be taken to 
be A^FP = BF X 156,097 = 0.0022, the product of the blend frequency for KOI-72.01 and 



the total number of Kepler targets (Borucki et al. 2010b). 



If we were to accept a confidence level of 3cr (99.73%) as sufficient for validation of 
a transiting planet candidate (corresponding to FAR = 2.7 x 10^'^), then the minimum 
number of super-earth-size planets required in order to be able to claim this level of 
confidence happens to be Ap = 1, according to the expression above. This value is so 
small that it gives us high confidence that the KOI-72.01 signal is not a false positive, but 
instead corresponds to a bona-fide super-earth-size planet. Another way to view this is 
that the expected number of background stars capable of producing the observed signal is 
~ 8 X 10^^ (see Table [sj bottom of column 5), so the probability that what we are seeing 
is a background blend is ~ 8 x 10^^ x Pgp/Pse, where Pop and Pse are the a priori 
probabilities that giant planets and super-earths occur with short-period orbits. Thus, 
claiming a FAR corresponding to a 3-cr detection only requires that the aforementioned 
probability ratio does not exceed 300. 

It is worth noting that the arguments above allow us to validate KOI-72.01 
independently of the detection of the refiex motion of the star (RV). This is possible 
in this case because we are able to rule out, using BLENDER, all false positives in which 
the background eclipsing binary is composed of two stars, which are predicted to induce 
ellipsoidal variation at a level that is not present in the photometry. Were this not the 
case, the blend frequency BF would have come out considerably larger, requiring in turn a 
significantly higher value for Np. 
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5.2. The KOI-72.02 Signal 

A similar BLENDER analysis was performed for KOI-72.02. For circular orbits, 
hierarchical triple systems are ruled out for the same reasons as in KOl-72.01. When 
considering the case of background eclipsing binaries (stellar tertiaries) with circular orbits, 
we find a range of blend scenarios that provide acceptable fits to the light curve, with 
secondaries of similar spectral type as the primary. Unlike the situation for KOI- 72. 01, 
the longer period of KOI-72.02 (45.3 days) leads to negligible ellipsoidal variation for the 
contaminating binaries, and this does not allow these types of blends to be excluded based 
on the quality of the fit, as we were able to do before. Chance alignment scenarios in which 
the tertiaries are planets rather than stars also lead to viable blends that can be up to 
4 mag fainter than the target. Those that are brighter than AKp = 2 can be ruled out 
because they would have produced a spectroscopic signature, but fainter ones cannot be 
ruled out any other way if they are angularly close enough to the target to be unresolved 
by our imaging observations. 

However, the longer period of this signal does not justify the assumption of a circular 
orbit. Allowing the orbit of a contaminating binary to have arbitrary eccentricity and also 
arbitrary orientation (longitude of periastron) can significantly increase the range of blends 
that provide good matches to the Kepler photometry, both for the chance alignment case 
and for physically associated triples. This is because the orbital speed in an eccentric orbit 
can be considerably larger or smaller than in a circular orbit, allowing for blends involving 
smaller or larger secondaries than would otherwise be permitted while still matching the 



observed duration of the transits, as described by Torres et al. (2010). It also increases 
the complexity of the problem, as the space of parameters to be explored is much larger. 
Additionally, because the secondaries are now not necessarily of the same spectral type as 
the primary star, attention must be paid to the resulting color of these blends, which could 
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be different from the measured colors as reported in the Kepler Input Catalog (Latham et 



al. 2005), in which case the blend would be excluded. Because these complications require 



significantly more effort to address, we are unable to provide sufficient evidence for the 
planetary nature of the KOI-72.02 signal at this time based on BLENDER considerations 
alone, and we defer such a study to a forthcoming publication. We note, however, that 
the transit duration and period of the two transit signals identified in the light curve give 
the same stellar density (1.142 ± 0.092 g cm'^ for KOI-72.01 versus 1.147 ± 0.096 g cm'^ 
for KOI-72.02) as derived from the transit properties using the analytic expression given 



in Equation 9 of Seager & Mallen-Ornelas (2003). This is a rare coincidence for a blend 
configuration. 



6. Stellar Characteristics 
6.1. Spectroscopic Parameters 

We carried out an LTE spectroscopic analysis using the spectral synthesis 



package SME (Valenti & Piskunov 1996 Valenti & Fischer 2005) applied to a high 



resolution template spectrum from Keck-HIRES of Kepler-10 to derive an effective 
temperature, T^s = 5705 ± 150 K, surface gravity, logg = 4.54 ± 0.10 (cgs), metalhcity, 
[Fe/H] = —0.15 ± 0.03, vsini = 0.5 ± 0.5 kms~^, and the associated error distribution 
for each of them. To refine the true parameters of the star, we took a novel path to 
constrain its surface gravity. The above effective temperature was used to constrain the 



fundamental stellar parameters derived via asteroseismic analysis (see Section 6.2). The 
asteroseismology analysis gave 4.341 ± 0.012 which is 0.2 dex lower than the SME value. 
The asteroseismology value is likely superior because of the high sensitivity of the acoustic 
periods to stellar radius. Still, the asteroseismology result depended on adopting the value 
of Tefj from SME. We recomputed the SME analysis by freezing (adopting) the seismology 
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value for log (7. This iteration yielded values of T^s = 5627 ± 44 K, [Fe/H] = —0.15 ± 0.04, 
and rotational vsini = 0.5 ± 0.5kms~^. The revised effective temperature was then put 
back into the asteroseismology calculation to further constrain the stellar radius and gravity. 
This iterative process converged quickly, as the log (7 from seismology yielded an SME value 
for Tgff that was only slightly different from the original unconstrained determination. 



We also measured the Ca II H&K emission (Isaacson & Fischer 2010), yielding a Mt. 
Wilson S value, 5=0.180 and logi?^^ = -4.89. Thus Kepler-10 is a magnetically inactive 
star, consistent with its low rotational rate, vsmi= 0.5 kms~^. Thus Kepler-10 appears 
to be an old (age greater than 5 Gyr) slowly rotating inactive star, slightly above the 
main sequence. This is consistent with the age derived from the asteroseismology analysis 
(Section [6^. 



6.2. Asteroseismology and the Fundamental Stellar Properties 

With a magnitude in the Kepler bandpass of Kp= 10.96, Kepler-10 presented itself 
as a promising case for asteroseismic characterization and was, consequently, placed on 
the short-cadence target list before it was even identified as a planet candidate. Figure 



TOfi illustrates the power density spectrum of the short-cadence light curve. It shows a 
clear enhancement of power, as expected for solar-like oscillations, around a frequency of 
2500 /iHz. In the spectrum one can identify sequences of approximately uniformly spaced 
peaks. This is in accordance with the asymptotic behavior of high-order acoustic modes, 
according to which the cyclic frequencies approximately satisfy 

Uni^^Mn + l/2 + e)-l{l + l)Do (1) 



( Vandakurov||1967 Tassoul||1980 ), where n is the radial order and / is the spherical-harmonic 



degree of the mode. The large frequency separation Vni — ^n-ii — Az/q is essentially given 
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by the inverse sound travel time across a stellar diameter; it is closely related to the 

1 /2 

mean stellar density p^, approximately satisfying Auq oc pj . For main-sequence stars 
Dq, giving rise to the small frequency separations Vni — ^n-i«+2, is largely determined by 
the variation of sound speed in the core of the star and hence provides a measure of the 
evolutionary state of the star. Finally, e is determined by conditions near the surface of 
the starH (For details on the diagnostic potential of solar-like oscillations, see, for example 



Christensen-Dalsgaard 2004). Photometric observations such as those carried out by Kepler 



are essentially restricted to degrees / < 2. 

The analysis of the observed frequency spectrum largely followed the procedures used 



by Christensen-Dalsgaard et al. (2010), with a pipeline developed for analysis of the Kepler 



p-mode data (Christensen-Dalsgaard et al. 2008 Huber et al. 2009). The first step was to 



carry out a correlation analysis to determine the large frequency separation Az/q, leading 
to Ai^o = 118.2 ± 0.2 yuHz. Using the scaling with the mean density and the corresponding 
values for the Sun, this yields a first estimate of = 1.080 ± 0.006 gcm"^. It should be 
noted, however, that this estimate does not take into account detailed differences between 
the structure of the star and the Sun. 

The next step in the analysis was to identify individual modes in the frequency 
spectrum. The detailed structure of the spectrum is illustrated in Figure [To| 3 which shows 
the folded spectrum, i.e., the sum of the power as a function of the frequency modulo Az/q- 
As indicated, in accordance with Eq. ([T]) there is clearly a pair of closely spaced peaks, 
corresponding to / = 2 and 0, as well as a single peak for / = 1. Given this identification, 
we were able to determine the individual frequencies of 19 modes. These are illustrated in 
an echelle diagram (cf. |Grec et 'al. 1983, see below) in Fig. 11 



^Owing to the small value of vsmi we do not have to consider rotational effects on the 
frequencies. 
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To determine the stellar properties, we fitted the observed frequencies to a grid of 
models. These were computed using the Aarhus Stellar evolution and pulsation codes 
(Christensen-Dalsgaard 2008a|[b ). See also Christensen-Dalsgaard et al. (2010). Diffusion 



and settling of helium and heavy elements were neglected. The grid consisted of models of 
mass between 0.8 and 1.1 Mq in steps of 0.02 Mq and mixing length oml = 1-5, 1.8 and 2.1. 
The composition was characterized by heavy-element abundances Z = 0.011,0.0127 and 
0.0144, corresponding to the observed [Fe/H] = —0.15 ± 0.06 (see above); this assumes a 
solar surface ratio between the heavy-element and hydrogen abundances {Zs/Xs)q = 0.0245 



(Grevesse & Noels 1993), and with the initial hydrogen abundance Xq related to Z, from 



galactic chemical evolution, through Xq = 0.7679 — 3Zq. 

For each model in the grid we calculated an evolution track extending well beyond the 
end of central hydrogen burning and, for the relevant models along the track, we computed 
adiabatic frequencies for modes of degree / = — 2. The match between the observed and 
computed frequencies, 1^1°!^^^ and was characterized by 

2 



1 / ..(mod) 1 



= - 

" N . 

nl 



where is the number of observed frequencies and ai, is the standard error in the observed 
frequencies, which we estimated as 1 /iHz. In the fit we also considered the observed 
effective temperature T^s (see above), using as combined measure of the goodness of fit 

(ry-i(obs) ry-,(mod) \ 

where T^^^^ and T^^°^^ are the observed and model values and a(Tes) is the standard error 
in Teff. 

Along each evolution track we determined that model which minimized x^? among the 
discrete timesteps in the evolution sequence. The final best model corresponding to the 
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track was obtained by further minimizing to determine Xmm for that track, interpolating 
Teff hnearly between timesteps and scahng the frequencies according to . Given the 
resuhing model values for all evolution tracks, we determined the final estimates of the 
stellar properties as an average over all tracks, weighted by XmL- 

A preliminary application of this procedure, using the value Teg = 5705 ± 150 K 
as described in the preceding section, resulted in an average surface gravity of 
logf? = 4.341 ± 0.012, substantially different from the spectroscopically determined value 
of 4.54 ± 0.10. In view of the known problems with the spectroscopic determination of 
logfyf we repeated the spectroscopic analysis, fixing \o^g as 4.341 ± 0.012 in accordance 
with the asteroseismic inference. This resulted in Tefj = 5627 ± 44 K. We then repeated the 
asteroseismic fits, conservatively increasing the uncertainty in Teg to 60 K. This resulted in 
the final estimates of the stellar parameters presented in Table [7j Mi, = 0.895 ± O.O6OM0, 
i?^ = 1.056 ± O.O21i?0, = 1.068 ± 0.008 g cm'^, and age = 11.9 ± 4.5 Gyr. The distance 
of 173 ± 27 pc is computed using the (/-band SDSS apparent magnitude in the Kepler Input 



Catalog and bolometric corrections interpolated form the tables of Girardi et al. (2008). 



To illustrate the quality of the fit. Figure [TT] shows the observed frequencies and 
the frequencies for the best-fitting model in an echelle diagram, reflecting the asymptotic 
structure of the spectrum described by Eq. ([T]). In accordance with this equation, the 
spectrum has been divided into segments of length Az/q which have been stacked. Formally 
this corresponds to reducing the frequencies modulo Auq. The asymptotic behavior is 
reflected in the nearly vertical columns of points, corresponding to the different values of 
the degree. Also, it is clear that the model provides a reasonable, although far from perfect, 
fit to the observations. E] 



''It may be noted that, unlike the simple expression (1), the small separation between 
modes with / = and 2 depends on frequency. 
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Planet Characteristics 



The physical and orbital properties of both transit signatures are derived by 
simultaneously fitting Kepler photometry and Keck RVs and by adopting the mean-stellar 
density of the host star as determined by asteroseismology. Our system model uses the 



analytic formalization of Mandel & Agol (2002 ) to fit photometric observations of the transit. 



We use the fourth-order non-linear parameterization of limb-darkening also described by 



Mandel fc Agol (2002) with coefficients (ci = 1.086, C2 = -1.366, C3 = 1.823, C4 = -0.672) 



calculated by Prsa (2010) for the Kepler bandpass. We account for variability phased to 



the orbital period by including the effects of reflected and emitted light from the planet, 
ellipsoidal variations due to tidal distortions of the host star and Doppler boosting due 
to motion of the star around the center of mass. For reflected/emitted light, we assume 
that the phased lightcurve is reproduced by a Lambertian reflector scaled by the geometric 



albedo. Ellipsoidal variations are modeled by using the prescription of Pfahl et al. (2008). 



Doppler boosting uses the methodology outlined in van Kerkwijk et al. (2010). We model 



the occultation by computing the fraction of the planet occulted by the star as a function 
of the star-planet projected distance. We assume that the planet is a uniformly illuminated 
disk during occultation. 

Our model parameters are the stellar mass and radius (M^, i?^), the planet mass and 
radius (Mp, -Rp), the orbital inclination (i), eccentricity (ecosw, esinw), the geometric 
albedo (^4^), and the RV amplitude and zero point {K, Vq). Model flts to the Kepler-lOb 
light curve yield an eccentricity that is consistent with zero {ecosw = 0.02 ± 0.10; 
esinw = —0.13 ± 0.20) which is consistent with our expectations for tidal circularization 



(Mazeh 2008). Given the large orbital separation of the outer planet candidate, we can 



not assume its orbit to be circular based on tidal circularization. The small predicted RV 
amplitude prevents a measurement of the eccentricity from existing RV observations. The 
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duration of the transit for the outer planet candidate is consistent with a circular orbit, but 
the resulting upper limit is still significant (e ~ 0.2). For the remainder of our discussion, 
the models are constrained to zero eccentricity for both Kepler-lOb and KOI-72.02. The 
radial velocity variations are modeled by assuming non-interacting (Keplerian) orbits. If 
KOI-72.02 were to be confirmed, then the relative inclination between the two orbits is 
likely less than 20°, as larger relative inclinations would require a fortuitous alignment of 



the orbital nodes for both planets to transit (Ragozzine & Holman 2010). 



We initially fit our observations by fixing and i?^ to their asteroseismic values 



see Section 6.2). Model parameters are found by chi-squared minimization using a 



Levenberg-Marquardt prescription. We then use the best-fit values to seed a Markov 



Chain Monte Carlo (MCMC) parameter search (Ford 2005) to fit all model parameters. 
We adopt the uncertainty on the asteroseismic determined mean-stellar density as a prior 
of the stellar mass and radius. A Gibbs sampler is used with widths initially defined by 



uncertainties derived from the diagonals of the constructed co- variance matrix (Ford 2006). 
Our Markov-chain contains 194,066 elements. We list the median of the distribution for 
each model parameter in Table [T] and the corresponding ±68.3% credible intervals (akin to 
a Icr confidence interval) centered on the median. The resulting properties of Kepler-lOb 
are as follows: Mp = 4.56;}:29 Rp='^-4^^Q-omlR(S, Pp = S-S^lg g cm-^, and a surface 
gravity (3.35^o dex) that is just 2.3 times that of Earth. Our knowledge of the planet is 
only as good as our knowledge of the parent star. Here, the planet radius is determined 
with a precision of just over 2% - comparable to the precision of the stellar radius derived 
from asteroseismology. The precision of the planet mass, however, is driven by the low 
semi-amplitude of the RVs (3.3!^° o nis~^) relative to the internal errors (1.5-2 ms~^). 

The absence of a statistically significant orbital signature in the RV data at the 
45-day period of KOI-72.02 translates to an upper limit for the mass under the planet 
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interpretation. The best fit to the RVs, constrained by the photometric period and phase, 
yields a shghtly negative semi-amphtude and, consequently, a mathematically valid but 
physically unrealistic negative planet mass: Mp = -2.80^^1^ Me- Of more relevance is the 
distribution of the masses returned by the MCMC calculations. The upper mass limit is 
taken to be three times the 68.3% credible interval (G.SM^), or 20 M^. This is the upper 
limit reported in Table [7| 

The phase-folded light curves together with the model fit are shown in the lower two 



panels of Figure 12 The best fit for Kepler-lOb requires a phase modulation with an 



amplitude of 7.6 ± 2.0 parts-per-million (ppm) and a 5.8 ± 2.5 ppm occultation, both of 



which are shown in the scaled and phase-shifted light curve in Figure 13 The modeled RV 
variations for Kepler-lOb are shown in Figure |6] as a function of orbital phase together with 
the 40 Doppler measurements (small circles) and the same averaged over 0.1-phase bins 
(large circles). 

The parameters derived from the MCMC analysis are listed in Table [7] and discussed in 
Section [9j Parameters such as mass, radius, and density for KOI-72.02 are included in the 
Table for completeness. However, the reader is cautioned against over-interpretation since 
the KOI-72.02 transits have not yet been confirmed to arise from a planetary companion. 



8. Statements about transit timing variations 

To measure the transit times, we generate a template transit shape based on folding the 
light curve with the given linear ephemeris. For each transit, we estimate the transit time by 
performing a local minimization, varying the transit mid-time, the light curve normalization 
and the light curve slope (outside of transit), but holding the remaining parameters 
(planet-star radius ratio, transit duration, impact parameter and limb darkening) fixed. 
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We iterate to improve the template transit shape used for measuring the transit times. 



The results are shown in Figure 14 The true uncertainties may be larger than the formal 
uncertainties (indicated by error bars), particularly for Kepler-lOb. Regardless, we do 
not detect statistically significant transit timing variations for Kepler- 10b at levels above 
0.01 days and for KOI-72.02 at levels above 0.003 days. Given the masses and periods we 
measured for both planet candidates, we predicted that TTVs would not be detectable, as 
indeed is the case. 

The precision with which the transit times of Kepler-lOb can be measured is lessened 
by the small size of the planet. Coupled with this loss in timing resolution, the period 
of the planet is also quite short at less than a day. Thus, the ratio of the period to the 
timing precision (the signal-to-noise ratio) is only P/a ~ (75000s)/(500s) = 150. For 
comparison, a Jupiter-size planet (~ 1% transit depth) in a 1-week orbit would produce 
a ratio about 100 times larger. The transit timing variations (TTV) from a nonresonant 
perturbing planet are essentially independent of the mass of the transiting planet. To 
summarize, additional non-transiting planets orbiting Kepler-10 are not well constrained 
by the measured transit times of Kepler-lOb. Here our RV measurements, which have 
very good precision of approximately 2.6 ms~^, are better suited to detect additional, 
nontransiting planets in this system. 

The small mass of this planet does provide good sensitivity to resonant perturbing 
planets where the TTV signal scales with the ratio of planet masses (see equation (33) of 



Agol et al. (2005)). Figure 15 shows the maximum allowed mass of a perturbing planet in 
an orbit near Kepler-lOb including both the transit time and RV data. Near mean-motion 
resonance, small planets with masses below that of Mars can be excluded. Away from 
resonance, the RV data constrain the presence of additional planets to be less than a few 
Earth masses. 
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If KOI- 72. 02 is, indeed, a longer-period planet orbiting Kepler-10, a TTV signal due to 
the interaction with Kepler-lOb must exist at some level and it behooves us to consider the 
expected magnitude of such a signal and whether or not it would be detectable in the data 
at hand. The hierarchical architecture limits the mechanisms that are capable of inducing a 
detectable signal. Neither changes in the light travel time caused by the displacement of the 
star due to a hypothetical outer planet, nor an evolving tidal field caused by an eccentric 
outer planet (see section 4 of Agol et al. ( 2005[ )), would be detectable in this system given 
the measured timing uncertainties. This statement is true even when one accounts for the 
factor of a/ Ntrans (whcrc Ntrans IS the number of observed transits) statistical improvement 
in the sensitivity. A scenario with an outer planet having a large eccentricity would only 
produce a detectable signal if its mass is comparable to Jupiter and its eccentricity is 



greater than 0.95 (see equation (25) of Agol et al. (2005)). Consequently, the absence of a 
TTV signal in the data does not rule out the planetary interpretation of KOI- 72. 02, given 
the upper mass limit of 20 (Section [T]). 



9. Discussion 



9.1. Composition of Kepler-lOb 



Kepler-lOb is a high-density rocky planet. This conclusion is based on the comparison 
of its radius and mass, measured within 1-sigma, with theoretical calculations of interior 



structure, as illustrated in Figures 16 and 17 The conclusion accounts for known 
uncertainties in the theory. 

Figure [16] is the mass-radius diagram for small planets in units of Earth mass and 



radius (after Zeng & Sasselov (2010)). Earth lies on the left edge, between two models (solid 
lines) of earth-like composition with a range of Fe/Si ratios. Uranus and Neptune are out of 
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range, above the upper right portion of the diagram. Kepler-lOb (bottom 1-sigma elhpse) 
and KOI-72.02 (1-sigma band with upper hmit in mass at ^ 25 Me) are shown with models 
and with two other planets known in this size range: CoRoT-7b and GJ1214b. In Figure 16 



the upper 1-sigma ellipses are for GJ 1214b (Charbonneau et al. 2009) and CoRoT-7b 



(Queloz et al. (2009) for mass; Bruntt et al. (2010) for radius). The uncertainty in the 
mass of CoRoT-7b is dependent on the methodology used to correct for the stellar activity 



signal that dominates the Doppler data. For example, Queloz et al. (2009) compute a mass 



of 4.8 ± 0.8 (the value used to construct the error ellipse for CoRoT-7b in Figure 16) 



while Ferraz-Mello et al. (2010) obtain 8.5 ± 1.5 M®. Hatzes et al. (2010) compute a mass 



of 6.9 ± 1.5 Me, while Pont, Aigrain, & Zucker (2010) obtain 2.3 ± 1.8 M®. The mass of 



CoRoT-7b, provided by Pont, Aigrain, & Zucker (2010), is shown by straight solid lines in 



Figure 16 in order to illustrate the range of mass values present in the literature. 



Theoretical calculations in Figure 16 are shown as curves, which are (from top to 



bottom): 10% by mass H/He envelope with typical ice giant interior similar to Uranus 
and Neptune (short-dashed line); theoretical pure water object (dot-dashed line); 50% 
water planets with 34% silicate mantle and 16% Fe core (thick long-dashed line), or with 
a low Fe/Si ratio of 44% mantle and 6% Fe core (thin long-dashed line; and earth-like 
composition with the same Fe/Si ratios (thick and thin solid lines). These models are 
from the grid by Zeng &: Sasselov| (2010), based on Valencia, O'Connell, & Sasselov (2007) 



with many updates, e.g. new water EOS by French et al. (2009). The Fe/Si ratios chosen 



in this grid correspond to the range observed in the Solar System given a stellar Fe/Si 



range in the solar neighborhood (Grasset, Schneider, & Sotin 2009), where the low Fe/Si 
ratio corresponds to the lowest values measured in small bodies in the outer Solar System, 
namely Ganymede. Note that the hypothetical 100% water and H/He envelope model 
curves are shown for illustration only; given the extremely high equilibrium temperature of 
Kepler- 10b an extended hot atmosphere of such volatiles has to be accounted for separately. 
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but Kepler-lOb is too dense for these scenarios. 

The interior structure of GJ 1214b has been modeled as an H/He/H20 planet with a 



rocky core (Nettelmann et al. 2010). Indeed, it lies between the models of ice giants similar 



to Neptune and Uranus and the models of a 50% water planet. And while the Queloz et 



al. (2009) mass and radius point to a rocky composition, the lower mass of Pont, Aigrain, 



& Zucker (2010) marginally favors a water/ice composition. The properties of Kepler-lOb 



together with their uncertainties clearly indicate a high-density, rocky planet. 



The dotted curve at the bottom of Figure [16] is the envelope corresponding to a 
maximum Fe core fraction expected from simulations of mantle stripping by giant impacts 



(Marcus et al. 2010). These simulations and planet formation scenarios indicate that pure 



Fe core objects cannot form in the mass range we consider. At the mass of Kepler-lOb 
these simulations predict that the maximum attainable Fe core is about 75% by mass. For 
comparison, Mercury has a ~70% Fe core. 

Theoretical calculations of the interiors of high-density planets suffer from a number 
of uncertainties in equations of states, high-pressure phases of materials inaccessible to lab 
experiments, cooling and differentiation histories, etc. There is consensus in the literature 



(Valencia, O'Connell, & Sasselov 2006, 2007 Fortney, Marley & Barnes 2007; Seager et 



al. 2007 Grasset, Schneider, & Sotin 2009) regarding the general results. However, the 



theoretical uncertainties are compounded by a degeneracy between only 2 observables 
(radius and mass), and 3 or 4 distinct types of bulk materials: Fe core, mantle (silicates. 



etc.), water, and hydrogen/helium gas (see Valencia, O'Connell, & Sasselov (2007); [Rogers 



fc Seager (2010) for details). 



At very high densities the radius of solid planets is constrained by the lack of bulk 
materials with density and compression properties above that of Fe and Fe alloys (see 



lower-envelope dotted curve in Figure 16). Planets with high densities close to that envelope 
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are likely composed predominantly of silicates and Fe; the degeneracy is lifted as indicated 
by the curvature of the iso-radius curves when they approach asymptotically the "dry" 
right-hand side of the ternary diagram in Figure IT7| 



Ternary diagrams were introduced to studies of solid planet structure by Valencia 



O'Connell, & Sasselov (2007). The three axes of the ternary diagram represent the core, 
mantle, and ice fraction of a planet of a given mass (in this case, the mass derived for 
Kepler-lOb). The axes are read by following lines parallel to the edges so that the three 
mass fractions sum to unity. Although there are 6 radiants emerging from each point in the 
grid, only three will result in mass fractions that add to unity. Different combinations of 
core, mantle, and ice fraction yield a planet radius that can be computed with theoretical 
models constrained by a given mass. Radius is not unique to a specific combination of mass 
fractions. The solid line in Figure [17] is an iso-radius. It shows the possible combinations 
that all yield the derived radius of Kepler-lOb. The dotted lines provide bounds on the 
domain captured by the one-sigma errors in planet radius. A ternary diagram represents a 
cross-section at a given mass. However, mass-dependent uncertainties are included in the 
error bars represented by the dotted line. The iso-radius together with the one-sigma error 
good representation of the possible bulk compositions for a given planet. 



The degeneracy in bulk composition discussed above is best illustrated in Figure [T7 
as a planet defined by a mass and radius is represented by a curve (solid) that can span a 
range of iron:mantle:water fractions. This degeneracy is practically lifted only for very high 
density planets near the "dry" right-hand side, close to the 100% iron core vertex, because: 
(1) the iso-radius curves bend as they approach the 0% water level, and (2) the H2O and 
OH molecules are able to be incorporated inside high-pressure silicate phases without 
changing their EOS much. For Kepler-lOb, the combination of very high density and 
constraints from mantle-stripping simulations (dash-dot line), restrict its bulk compositions 
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to dry rocky iron-core-dominated interiors similar to Mercury. On the other hand, if the 
observational derivation of the mass of Kepler- 10b were significantly overestimated (by 
more than 1-sigma), the planet could contain significant amount of water in its interior. 

In conclusion, within l-to-2-sigma in its derived radius and mass, the planet Kepler-lOb 
is a dry rocky planet with high Fe content. Its high density does not violate the prediction 



for maximum mantle stripping during planet formation (Marcus et al. 2010) 



9.2. The Phase Curve of Kepler-lOb 

The phase curve amplitude of Kepler-lOb is 7.6 ± 2.0 ppm (See Table [?]). If due to 
scattered light alone, this corresponds to a Kepler bandpass effective geometric albedo of 
0.68. The occultation depth is 5.8 ± 2.5 ppm which corresponds to an effective geometric 
albedo of 0.61 ± 0.17. This is an unusually high albedo. The only known solar system 
bodies that are so bright are Venus (due to photochemically induced hazes) and Saturns icy 
moon Enceladus (coated with fresh ice). Kepler-lOb is likely too hot for any hazes and is 
certainly too hot for any ice. Another possibility is that Kepler-lOb has silicate clouds, but 
they would have to be completely covering the planet's day side and have a large particle 



size in order to provide the required reflectivity (see, for example, Seager, Whitney, & 



Sasselov (2000)). 



We prefer an interpretation that the phase curve is dominated by a thermal radiation 
change from the planet's day to night side. This case would be reminiscent of hot Jupiters, 
which have a hotter day side than night side. In the case of the Kepler bandpass, the 
temperature difference need not be too extreme, as long as the planet is hot enough on the 
day side (as is the case for Kepler-lOb with a equilibrium temperature of 1833 K). This is 
because the Kepler bandpass is in the optical, and so the contribution of thermal radiation 
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drops off rapidly with decreasing temperature (i.e. day-to-night side). In other words, a 
passband on the Wien side of the black body curve will capture significantly different ffuxes 
with just a small change in temperature. An example of this is the phase curve of HAT-P-7 



Borucki et al. 


(2009 


); 


Welsh et al. 


(2010 



An intriguing possibility is that Kepler-lOb has no atmosphere at all, having lost 
it over time due to atmospheric erosion. One test is to search for a hotspot directly at 
the substellar point, since no atmospheric winds would be available to move the hotspot 



off center (Seager & Deming 2009). Further Kepler data, in particular a confirmation 
and robust measurement of the occultation depth, of Kepler-lOb will help with any 
interpretation. 

Asteroseismic measurements of the host star Kepler-10 indicate that the planetary 
system is very old, allowing for more than 11 Gyr of evaporation and mass loss from the 
surface of Kepler-lOb. Assuming, as is common, that the planet arrived at its present orbit 
within about 100 Myr of formation, the relevant evaporation time scales are those of water 
steam or a silicate surface. Any H/He envelope would have evaporated too quickly (at ~ 
10^^ g s^^) to affect the further evolution of the planet, while the evaporation of the silicate 



mantle is too slow to remove more than 50% of the planet's original mass (Valencia et al. 



2009). Therefore, it is difficult to establish if Kepler-lOb is the remnant core of a water 
planet or an ice giant planet, though it is probable that it must have lost a significant 
fraction of its mass. 



10. Summary 



NASA's Kepler Mission collects transit photometry from a spaceborne Schmidt camera 
to detect and characterize extrasolar planets with the goal of determining the frequency of 
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earth-size planets in or near the habitable zone of Sun-like stars. Now in its second year 
of operation, the mission has reached an important milestone toward meeting that goal, 
namely the discovery, reported herein, of its first rocky planet, Kepler-lOb. This planet was 
identified via transit photometry. The analyses described here are based on ~ 8 months of 
Kepler 29.4-minute cadence data acquired between 02 May 2009 and 09 January 2010. The 
target was also observed at a higher 1-minute cadence from 21 July 2009 to 19 August 2009 
and 18 September 2009 and 09 January 2010. 

Two distinct sets of transit events were detected in the lightcurve of Kepler-10 
constructed from ~ 8 months of Kepler photometry: 1) a 152 ± 4 ppm dimming 
lasting 1.811 ± 0.024 hours with transit ephemeris of T[BJD] = 2454964.573751^;^ + 
N * 0.8S7A95to!S)t days and 2) a longer-period event described by a 376 ± 9 ppm 
dimming lasting 6.86 ± 0.07 hours and an ephemeris T[BJD] = 2454971.676112:0^23 + 
N * 45.29485lo;ooo76 days. Statistical tests on the photometric and individual pixel flux 
time series of Kepler-10 estabhshed the viabihty of the planet candidates. For example, 
comparison of the flux-weighted photocenter during transit and outside of transit revealed 
no deviation consistent with an eclipse event associated with one of the nearby stars 
that might be diluting the lightcurve of Kepler-10. Clean statistics triggered a battery of 
ground-based follow-up observations. 

High resolution reconnaissance spectroscopy was used to verify the effective temperature 
and surface gravity as well as rule out obvious eclipsing binaries masquerading as planets 
by way of moderate-precision RVs. High spatial resolution imaging (AO and Speckle) 
was acquired to identify faint, nearby stars that should be considered in the photocenter 
analysis as potential background eclipsing binaries. In the case of Kepler-10, no additional 
stars were identifled. Forty precision Dopplcr measurements were acquired with the Keck 
10-meter telescope between May and August 2009. These measurements conflrmed the 
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planetary nature of the short-period transit event. The photometric period was clearly 
seen in a periodigram of the velocities, and the variations are phased as expected given the 
transit epochs. With a semi- amplitude of just 3.3^5!o the Doppler measurements 

suggest a planetary mass for this companion. No significant signal was detected in the 
measurements at the photometric period of the outer candidate, KOI- 72. 02. 

Knowledge of the planet is only as good as our knowledge of the star it orbits. Matching 
the reconnaissance spectroscopy of Kepler-10 to a library of synthetic spectra yielded 
Teff = 5680 ± 91 K, log^ = 4.33 ± 0.16, \Fe/H] = -0.09 ± 0.04, and vsmi = 1.5 ± 0.5 
kms~^. Full spectral synthesis using HIRES echelle data without the iodine cell yielded 
Teff = 5705 ± 150 K, log^ = 4.54 ± 0.10, [Fe/H] = -0.15 ± 0.03, and vsini ^ 0.5 ± 0.5 
kms"-*^. Because the parent star is relatively bright (Kp = 10.96) and, hence, amenable 
to asteroseismic analysis, Kepler photometry was also collected at 1-minute cadence for 
~ 5 months from which we detected 19 distinct pulsation frequencies. Modeling of these 
frequencies resulted in precise knowledge of the fundamental stellar parameters. The 
process was iterated once in that the asteroseismic analysis yielded a surface gravity of 
log g — 4.341 ± 0.012 which was then fixed in the spectral synthesis analysis to yield 
an improved effective temperature, — 5627 ± 44 K, that was then fed back to the 
asteroseismic analysis. The result is that Kepler-10 is a relatively old (11.9 ± 4.5 Gyr) but 
otherwise Sun-like Main Sequence star with Teg = 5627 ± 44 K, = 0.895 ± 0.060 M©, 
and = 1.056 ±0.021 Rq. 

Physical models, constrained by the asteroseismology-derived stellar parameters, were 
simultaneously fit to the transit light curves and the precision Doppler measurements. 
Modehng produced tight constraints on the properties of Kepler-lOb: Mp — 4.56^^9 
M®, i?p=1.416lQ;g3gi?0, andpp = 8.8l|g g cm~^. Evaluation of these properties within 
a theoretical framework allowed us to draw conclusions about the planet's composition. 
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Within 1-2(7 of the derived mass and radius, Kepler-lOb is a dry, rocky planet with high Fe 
content. Its high density does not violate predictions for maximum mantle stripping during 
planet formation. 

The authors would like to thank Carly Chubak for computing the barycentric radial 
velocity of Kepler-lOand Frederic Pont for his careful reading and thoughtful comments, 
providing us with feedback in record time. JCD acknowledges support from The National 
Center for Atmospheric Research which is sponsored by the National Science Foundation. 
Funding for this Discovery mission is provided by NASA's Science Mission Directorate. 
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Table 1. Observed flux-weighted centroid shifts for KOI- 72. 01. 



Ql Q2 Q3 

AR 6.07 X 10-5 ± 2.42 x IQ-^ 7.36 x 10"^ ± 2.40 x 10'^ -4.91 x lO"'' ± 2.42 x 10"^ 
AC 6.96 X 10-5 ± 1.59 x IQ-^ -1.12 x lO'^ ± 1.25 x 10"^ -8.75 x IQ-^ ± 1.59 x lO'^ 
D 9.23 X 10-5 ± 1.99 X IQ-s 1.34 x 10-5 ± 9.25 x 10-5 8.77 x 10'^ ± 1.59 x 10-5 
D/a 4,64 ora 055 



Table 2. Observed flux-weighted centroid shifts for KOI-72.02. 



Ql 



Q2 



Q3 



AR 2.70 X IQ-"' ± 4.08 x IQ-'' 
AC 7.58 X 10-^5 ± 2.87 x 10"^ 
D 2.81 X 10-" ± 4.01 X 10-5 
D/a 7m 



-1.90 X IQ-'' ± 5.35 X IQ-"' 
1.91 X 10""' ±4.03 X 10^5 
2.69 X 10-* ±4.73 X 10-5 
5.68 



3.14 X 10-5 ± 3.05 X 10-5 
-3.84 X 10-5 ± 2.45 x 10-5 
4.96 X 10-5 ± 2.70 X 10-5 
1.83 



Table 3. PRF-fit centroid shifts for KOI-72.02 without saturated pixels. 



Ql 



Q2 



Q3 



AR 
AC 

D 

D/a 



3.05 X 10-" ± 7 X 10-5 
1.72 X 10-" ± 7 X 10-5 
3.50 X IQ-* ± 7 X 10-5 
4.9 



-1.82 X 10-" ± 7 X 10-5 
-1.84 X 10-" ± 7 X 10^5 
2.59 X 10-" ± 7 X 10-5 
3.7 



5.50 X 10-5 ±7 X 10-5 

4.39 X 10-5 ± 7 X 10-5 
7.05 X 10^5 _|_ 7 X 10-5 
1.0 
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Table 4. Modeled centroid shifts due to transits on the known stars in the Q3 aperture 
with depths that reproduce the observed Q3 depth. 



Object Modeled depth Modeled D D/a Object Modeled depth Modeled D D/a 



KOI-72.01 


1.67 X 10" 


4 


6.10 X 


10" 


-6 


0.384 


KOI-72.02 


5.52 X 10" 


4 


2.02 X 


10- 


5 


0.747 


11904143 


1.90 X 10" 


1 


5.69 X 


10" 


-4 


35.8 


11904143 


6.30 X 10" 


1 


1.89 X 


10" 


3 


69.7 


11904165 


6.31 X 10"^ 


2 


3.53 X 


10" 


-4 


22.2 


11904165 


2.09 X 10" 


1 


1.17 X 


10- 


3 


43.3 


11904167 


1.17 X 10" 


1 


4.43 X 


10" 


-4 


27.9 


11904167 


3.88 X 10" 


1 


1.47 X 


10- 


-3 


54.3 


11904169 


7.13 X 10" 


1 


6.73 X 


10" 


-4 


42.3 


11904169 


2.36 












11904171 


2.22 X 10" 


1 


6.94 X 


10- 


-4 


43.6 


11904171 


7.35 X 10- 


1 


2.30 X 


10- 


-3 


84.9 


11904145 


1.55 X 10" 


1 


4.04 X 


10- 


-4 


25.4 


11904145 


5.13 X 10- 


1 


1.34 X 


10- 


-3 


49.5 


11904150 


1.10 












11904150 


3.63 












11904152 


1.95 X 10" 


1 


5.46 X 


10- 


-4 


34.3 


11904152 


6.46 X 10- 


1 


1.81 X 


10- 


3 


66.8 


11904154 


3.08 X 10" 


1 


6.51 X 


10- 


-4 


40.9 


11904154 


1.02 












11904155 


1.56 












11904155 


5.15 












11904158 


3.09 X 10" 


2 


6.53 X 


10- 


-4 


41.1 


11904158 


1.02 X 10- 


1 


2.16 X 


10- 


-3 


80.0 


11904159 


1.16 X 10" 


1 


5.63 X 


10- 


-4 


35.4 


11904159 


3.82 X 10- 


1 


1.86 X 


10- 


3 


69.0 


11904160 


2.51 X 10" 


2 


2.04 X 


10- 


-4 


12.8 


11904160 


8.31 X 10- 


2 


6.74 X 


10- 


4 


24.9 


11904162 


1.06 X 10" 


1 


4.09 X 


10- 


-4 


25.7 


11904162 


35.1 X 10- 


1 


1.35 X 


10- 


3 


50.1 



Note. — Transits on some companions can be ruled out because they require depth > 1. 
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Distance Distance AJ J 
(FWHM) (arcsec) (mag) (mag) 



1 


0.075 


1.5 


11.4 


2 


0.150 


3.5 


13.4 


3 


0.225 


5.0 


14.9 


4 


0.300 


5.0 


14.9 


5 


0.375 


5.5 


15.4 


6 


0.450 


5.5 


15.4 


7 


0.525 


6.5 


16.4 


8 


0.600 


7.0 


16.9 


9 


0.675 


7.5 


17.4 


40 


3.000 


9.5 


19.4 



Table 5: Palomar AO source sensitivity as a function of distance from the primary target at 
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Table 6. Relative Radial Velocity Measurements of KOI- 72 . 



HJD RV ctrv 

(-2450000) (ms-^) (ms-^) 



5074.878 


4.06 


1.6 


5075.773 


6.83 


1.5 


5076.863 


0.40 


1.5 


5077.923 


-6.36 


1.5 


5078.922 


1.76 


1.7 


5079.973 


-2.41 


1.7 


5080.896 


5.34 


1.5 


5081.969 


-7.16 


1.5 


5082.848 


-8.50 


1.5 


5083.761 


3.27 


1.4 


5083.945 


-2.27 


1.6 


5084.878 


1.17 


1.4 


5106.890 


-0.14 


1.7 


5169.725 


1.58 


0.9 


5170.725 


-3.00 


1.1 


5172.756 


5.05 


1.6 


5173.721 


4.90 


1.1 


5312.047 


-2.68 


1.6 


5313.004 


-5.29 


1.5 


5314.005 


-9.82 


1.6 
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Table 6 — Continued 



HJD RV (Trv 

(-2450000) (ms-i) (ms-^) 



5317.998 


3.92 


1.7 


5318.121 


-0.70 


1.5 


5319.027 


-5.43 


1.6 


5320.063 


1.44 


1.5 


5321.007 


0.97 


1.6 


5321.969 


1.61 


1.5 


5343.050 


-4.01 


1.5 


5344.032 


-4.93 


1.5 


5344.973 


-4.08 


1.3 


5345.068 


-2.08 


1.7 


5350.973 


0.56 


1.6 


5351.988 


0.88 


1.4 


5373.814 


-0.19 


1.6 


5376.865 


-0.38 


1.4 


5379.902 


1.78 


1.7 


5403.898 


0.66 


1.4 


5407.013 


-3.80 


1.5 


5411.986 


-3.83 


1.4 


5412.805 


-4.36 


1.2 


5414.803 


-1.20 


1.5 
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Table 6 — Continued 



HJD RV (7RV 

(-2450000) (ms-^) (ms"^) 
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Table 7. Star and planet parameters for the Kepler- 10 system. 



Parameter Value Notes 



Transit and orbital parameters: Kepler-lOb 

Orbital period P (days) 0-837495tg;ggggg^ A 

Midtransit time E (BJD) 2454964.57375^0 ooo82 ^ 

Scaled semimajor axis a/R-^ 3.436^q Q92 ^ 

Scaled planet radius i?p/R* 0.012321°;°°°}^ A 

Impact parameter 6 0.339lo g7g A 

Orbital inclination i (deg) 84.4?olJ;g A 

Orbital semi-amplitude JC (ms-i) 3.3t?;^ B 

Orbital eccentricity e OB 

Center-of-mass velocity 7 (ms-i) -98.93 ± 0.02 B 

Transit and orbital parameters: KOI-72.02 

Orbital period P (days) 45.29485t!!;"™57g A 

Midtransit time E (HJD) 2454971.676llo;gg23 A 

Scaled semimajor axis a/R* 49.1^2'3 A 

Scaled planet radius Rp/R* 0.0193781o;ooom A 

Impact parameter b 0.299j]Q Qy| A 

Orbital inclination i (deg) 89.7?0l°;?| A 

Observed stellar parameters 

Eflfective temperature Tefi(K) 5627 ± 44 C 

Spectroscopic gravity logp(cgs) 4.35 ± 0.06 C 

Metallicity [Fe/H] -0.15 ±0.04 C 

Projected rotation v sin i(kms-i) 0.5 ±0.5 C 

Fundamental Stellar Properties 

Mass M*(M0) 0.895 ±0.060 D 

Radius _R*(_R0) 1.056 ±0.021 D 

Surface gravity logc/* (cgs) 4.341 ± 0.012 D 

Luminosity L* (Lq) 1.004 ± 0.059 D 

Absolute V magnitude My (mag) 4.746 ± 0.063 D 

Age (Gyr) 11.9 ±4.5 D 

Distance (pc) 173 ±27 D 

Planetary parameters: Kepler-lOb 

Mass Mp (Me) 4.56lJ;^^ A,B,C,D 
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Table 7 — Continued 



Parameter 


Value 


Notes 


Radius Rp {R^) 


1 416+""^^ 

l.'±lU_g Q3g 


A,B,C,D 


Density pp (gcm~^) 


».»_2 9 


A,B,C,D 


Surface gravity loggp (cgs) 


•^••^^-0.13 


A,B,C,D 


Orbital semimajor axis a (AU) 


U.Uib84_Q 0QQ34 


E 


Equilibrium temperature Teq (K) 


1833 


F 


Parameters for candidate: KOI-72.02 






Mass Mp (Me) 


< 20 


G 


Radius Rp (-Rgj) 


9 997+0.052 
_0.057 


A,D 


Orbital semimajor axis a (AU) 


U.^4U^_Q QQ53 


E 


Equilibrium temperature Teq (K) 


485 


F 



Note. — A: Based primarily on an analysis of the photometry, 
B: Based on a joint analysis of the photometry and radial velocities, 
C: Based on an analysis by D. Fischer of the Keck/HIRES template 
spectrum using SME i jValenti fc Piskunov|1996^ , 
D: Based on asteroseismology analysis E: Based on Newton's revised 
version of Kepler's Third Law and the results from D, 
F: Calculated assuming a Bond albedo of 0.1 and complete redistri- 
bution of heat for reradiation. G; Upper limit corresponding to three 
times the 68.3% credible interval from MCMG mass distribution. 
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Table 8. Blend frequency estimate for KOI-72.01 based on the frequencies of transiting 

giant planets. 



Kp range 




Stellar density 


pmax 


Stars 


Transiting Jupiters 


Transiting Neptunes 


(mag) 


(mag) 


per sq. deg 


( ) 


(x lO*") 


6-15 iJ®, /jup = 0.117o 


3.8-6 K®, /Nop = 0.07470 












(xlO-^) 


(xlO-'5) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


11.0-11.5 


0.5 












11.5-12.0 


1.0 












12.0-12.5 


1.5 












12.5-13.0 


2.0 












13.0-13.5 


2.5 


48 


0.12 


0.168 


0.0002 


0.0001 


13.5-14.0 


3.0 


87 


0.15 
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0.0005 


0.0003 


14.0-14.5 


3.5 
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0.18 
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0.0009 


0.0006 


14.5-15.0 
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0.20 
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0.0014 


0.0009 


15.0-15.5 


4.5 


189 


0.22 


2.217 


0.0024 


0.0016 


15.5-16.0 


5.0 


185 


0.25 


2.803 


0.0031 


0.0020 


16.0-16.5 


5.5 












16.5-17.0 


6.0 












Totals 




746 




7.766 


0.0085 


0.0055 



Blend frequency (BF) = (0.0085 + 0.0055) X lO"*' = 1.4 X 10" 



Note. — Column 1: magnitude bins; column 2: magnitude difference l\Kp compared to the primary, taken at the 
upper edge of each bin; column 3: density of stars from the Besangon models; column 4: angular separation at which 
stars in the corresponding magnitude bin would go undetected in our imaging observations; column 5: number of 
stars in a circle of radius pmax around KOI- 72 , restricted also by mass according to the contours in Figure [s] column 
6: number of expected transiting Jupiters as contaminants; column 7; number of expected transiting Neptunes as 
contaminants. 
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Fig. 1. — Raw (upper) and corrected (lower) flux time series produced by the Kepler pho- 
tometry pipeline (PA and PDC fluxes, respectively). Vertical lines denote the boundaries 

between quarters. Intra-quartcr flux values have been scaled by their median flux for display 
purpose only (to mitigate the discontinuities between quarters). 
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Fig. 2. — Speckle reconstructed R band image of Kepler-10. No other source is observed to 
within 1.8 arcseconds of the target to a depth of ~6 magnitudes. The cross pattern is an 
artifact of the reconstruction process. North/East are up/left. 
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Fig. 3. — J-band Palomar adaptive optics image of KOI-72 . The top image displays a 
12.5" X 12.5" field of view centered on the primary target. The bottom image displays a 
2" X 2" field of view centered on the primary target. The four-point pattern surrounding the 
central point spread function core is part of the adaptive optics point spread function. 
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Fig. 4. — Radial velocities derived from HIRES spectra collected in 2009 and 2010 are plotted 
against time. Error bars include not only the expected instrumental noise but also a 2 ms~^ 
jitter to account for variations intrinsic to the star. 
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Fig. 5. — A periodogram of the radial velocities shows a peak in the spectral density at the 
photometric period derived from the transits of Kepler-lOb. 
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Fig. 6. — Radial velocities versus phase derived from transit photometry of the short-period 
event. Both individual velocities are plotted (small circles) as well as averages over 0.1-phase 
bins (large circles). The dashed line shows the best-fit circular orbit solution for which there 
are only two free parameters - the amplitude K and the zero-point of the velocities. 
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Fig. 7. — Line bisector span measurements folded at the photometric period. 
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Fig. 8. — Map of the surface corresponding to a grid of blend models for KOI-72.01 
involving background eclipsing systems in which the tertiary is a (dark) planet, in a circular 
orbit around the secondary. Differential extinction is included (see text). The vertical 
axis represents a measure of the relative distance between the background binary and the 
primary star, which we parametrize here for convenience in terms of the difference in distance 
modulus. Contours are labeled with the Ax^ difference compared to the best planet model 
fit (expressed in units of the significance level of the difference, cr). Two dashed lines are 
also shown that correspond to equal magnitude difference (Aii'p) between the contaminating 
background binary and the primary star. Kinks in the contours are simply a result of the 
discreteness of the grid. 
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Fig. 9. — Sensitivity to faint companions near KOI-72 from our imaging observations. Any 
companions above the curves would be bright enough to be detected. J-band hmits are from 
AO observations at the Palomar 200-inch telescope, and R is from speckle observations using 
the WIYN 3.5 m telescope. The vertical dotted lines indicate the 3-cr confusion radius of 
1'.'17 for KOI- 72. 01 and 0'.'60 for KOI-72. 02, i.e., the maximum angular separation at which 
background eclipsing binaries would remain undetected in our centroid motion analysis. 
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Fig. 10. — Panel a) shows the power-density spectrum, with a 1.6 //Hz smoothing, for short- 
cadence observations during the second month of Q2 and all of Q3. Panel b) shows the 
spectrum folded at the large separation frequency of 118.2 //Hz where the relative power 
excess is the power minus the noise floor, the result of which is divided by the peak power. 
This allows identification of the peaks corresponding to modes of degree I — 0, 1 and 2, as 
indicated. 
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Fig. 11. — Echelle diagram (cf. 



Grec et al. 



1983 ) illustrating the observed frequencies (filled 



symbols) and the frequencies of one of the best-fitting model (open symbols), for modes of 
degree I = (circles), I = 1 (triangles) and I = 2 (squares). The model has a mass of 0.9 Mq, 
Z = 0.0144, auL = 1-8 and an age of 11.6 Gyr. 




Fig. 12. — The Kepler photometry and physical models are plotted as a function of both time 
(upper two panels) and phase (lower two panels). The transits of KOI- 72. 02 are highlighted 
by blue vertical lines in the topmost panel while a cutout (defined by the yellow box) is 
expanded to show the transits of Kepler-lOb highlighted by red vertical lines (second from 
top). The bottom two panels shows the phase-folded light curves centered on phase zero 
as defined by the central transit time. The modeled light curves are shown as colored lines 
(blue corresponding to KOI-72.02 and red corresponding to Kepler-lOb. Also shown is a 
phase cutout of the lightcurve and model (green) centered on phase= 0.5 where occulations 
would occur for a circular orbit. The relative intensity scale for phase= can be read off 
the y-axis on the left- hand- side of the plot while the relative intensity scale for phase= 0.5 
can be read off the y-axis on the right-hand-side of the plot. 



-77- 




Fig. 13. — The relative intensity scale of the phase-folded light curve of Kepler-lOb shown 
in the bottom panel of Figure 12 is expanded to show the phase modulation and marginal 
occultation required by the model fits. Colors have the same meaning as in Figure 12 
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Fig. 14. — The difference between eacfi best-fit transit time and a linear epfiemeris for 
Kepler-lOb (top) and KOI-72.02 (bottom). As expected, we do not detect any statistically 
significant transit timing variations. 
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Fig. 15. — Constraints (95% confidence level) on low-eccentricity secondary planets that are 
exterior to Kepler-lOb as a function of the period ratio using the measured transit times and 
RV data. The dotted curve shows the limit from a TTV analysis alone from equation (A7) in 



Agol et al. (2005 ). The dashed line is the expected sensitivity from 40 RV measurements with 



2.6 m/s precision calculated using equation (2) from Steffen & Agol (2005). The solid curve 
is the overall sensitivity from both RV and TTV measurements (summed in quadrature). 



The diamonds are calculations for MMR from equation (33) in Agol et al. (2005). The TTV 
sensitivity curve has been scaled down by \/269 to represent the improvement in sensitivity 
due to the number of observed transits. Finally, the horizontal dot-dashed and triple-dot- 
dashed lines correspond to the mass of the Earth and the mass of Mars respectively. 
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Fig. 16. — Kepler-lOb and KOI-72.02 on the Mass-Radius diagram. Kepler-lOb is shown with 
a 1-sigma elhpse at Rp = 1.416 Re] KOI-72.02 is shown with a 1-sigma band, constrained 
in mass below ^25Me- Two similar planets are shown for comparison: GJ1214b (1-sigma 
ellipse top) and CoRoT-7b (1-sigma ellipse below, and solid crossed lines for the alternative 
mass estimate as described in text). Theoretical models are shown as curves (from top to 
bottom): 10% by mass H/He envelope with typical ice giant interior similar to Uranus and 
Neptune (short-dashed line); theoretical pure water object (dot-dashed line); 50% water 
planets with 34% silicate mantle and 16% Fe core (thick long-dashed line), or with a low 
Fe/Si ratio of 44% mantle and 6% Fe core (thin long-dashed line); and earth-like composition 



with the same Fe/Si ratios (thick and thin solid lines). Models are from the grid of Zeng & 



Sasselov (2010). The dotted curve at the bottom corresponds to a maximum Fe core fraction 



expected from simulations of mantle stripping by giant impacts (Marcus et al. 2010). 
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Fig. 17. — The interior structure models for Kepler-lOb shown on a ternary diagram, il- 
lustrating the planet's solid nature. Three possible bulk materials can mix to determine a 
planet's radius at a fixed planet mass. The model planet radius appears as a curved line; 
this iso-radius curve for Kepler-lOb (solid line) lies very close to the "dry side" of the ternary 
diagram, indicating that most models at this high density exclude water as a significant bulk 
component of the planet's composition. The dash-dot line corresponds to a maximum Fe 



core fraction expected from simulations of mantle stripping by giant impacts (Marcus et al. 
(|20l0|) . 



